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1.0 SUMMARY 


This design note contains a description of the modifications necessary to 
give the Houston Operations Predictor/Estimator (HOPE) program the capability 
to solve for or consider vent forces. The HOPE program, with the above 
modifications, will be referred to in this paper as HOPE/VENT. A new 
HOPE version which will include HOPE/VENT as well as other recent program 
modifications is currently being prepared and will be released as HOPE/MDC 3.0. 

The formulation of HOPE/VENT relies heavily on the venting capability which 
already exists in version MDC 2.0 of HOPE. The user must input an attitude 
timeline as well as vent on-off times. The program will solve for or consider 
the components of a vent in body coordinates. This approach conforms to the 
current HOPE formulation as well as the approach taken in the Real Time 
Computer Complex (RTCC) logic. 

The verification of HOPE/VENT was naturally divided into two phases: solve 
for and consider capability. The solve for capability was tested using the 
existing venting logic in the dummy data mode. The consider capability was 
verified by comparing results with the H0PE/B8.02 IGS Burn Model. In the 
process of verifying the consider capability, numerous difficulties were 
•encountered with the HOPE trajectory integrator. For the proper functioning 
of the consider vent force model, the present study suggests a modification 
of the integrator to invoke a restart when encountering vent force 
discontinuities. 


Section 3.0 describes in detail the model Implemented in solving for vent 
forces, while Sections 3.1 and 4.2 detail the integrator problems 
encountered. Appendix A provides a summary derivation of the mathematical 
principles applicable to solve/consider methodology. 
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2.0 INTRODUCTION 


One of the many unique problems encountered in performing orbit determination 
processing for the Shuttle is the existence of unknown perturbative forces on 
the vehicle caused by venting various gases during orbital flight. Compli- 
cating the relatively simple effect of venting is the phenomenon of plume 
impingement. As gases are released, they interact with the various body 
surfaces of the vehicle, causing unknown accelerations. While the primary 
venting effect may be known beforehand, the secondary interaction with tail, 
wing and door surfaces is unknown. An additional effect is the existence of 
RCS uncoupled thrusting which is triggered by vent forces. 

The method described in this design note, called the method of correction 
vents, presumes the existence of a nominal vent timeline which must be 
corrected by solving for unknown vent perturbations. Thus, the nominal 
timeline may represent the primary effect of venting while the correction 
vents correspond to plume impingement effects. 

An active vent timeline may have forty (40) distinct entries in the HOPE vent 
table. Computer time limitations and questionable mathematical validity, 
however, suggest that solving explicitly for such a large number of vents 
would be impractical. Using the method of correction vents, instead, enables 
one to solve for a long term average vent error which may, in a least square 
error sense, realistically account for the effect of many error vents. 

In implementing the method of correction vents, an attempt has been made to 
keep program changes small. Rather than writing entirely new subroutines, 
use has been made of existing code whenever possible. Also, the unusual step 
has been taken to adapt little-used Lunar MASCON internal storage logic in 
the current HOPE for implementation of the new vent force capability. 


2 


The resulting loss In MASCON computational capability was considered a 
justifiable tradeoff for having an up-to-date Shuttle era program developed 
per the established venting modifications schedule. Moreover, plans have 
been made for reihcorporation of the MASCON capability (albeit without, 
simultaneous execution of the venting calculations), pending completion 
of the present vent modification performance studies. 
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3.0 DISCUSSION 

The approach being implemented can be best described as the method of cor- 
rection or error vents. It Is presumed that the user will input an attitude 
timeline and a table of nominal vents. The nominal vents are input in the 
body axis frame of the vehicle and should represent the user's best estimate 
of venting activity. The nominal vent timeHne is input through the logic 
which already exists in HOPE version MDC 2.0 (References 1 and 2). 

The vents to be solved for or considered, i.e., the correction vents, are viewed 
as updates or refinements to the nominal vent timeline. The correction vents 
are also input in the body axis frame and the steps necessary to do so are 
described in the user's section of this paper. Instead of solving for each 
nominal vent individually, the HOPE/VENT program solves for or considers an 
average vent error vector which can stretch over a fairly extended period of 
time. The number of such average error vents is limited; also, the user is 
required to input on-off times corresponding to each correction vent. The 
on-off times of the correction vents might represent special events in the 
vehicle attitude timeline, intervals of known extended vent activity, or 
other a priori vent information which the user might care to incorporate. 


The following example illustrates the HOPE/VENT program nomenclature just 



Vehicle Attitude 
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In the figure, the solid lines represent magnitudes of nominal vents, vj, v 2 , 
V3, and v 4 . The dashed lines represent correction vents, e^, e 2 . A point 
which should be stressed is that because of the computer implementation, 
start-stop times of correction vents must correspond to start and stop times 
of vents within the nominal table. This point is reiterated in the user's • 
section. In this example, the correction vents are put at the boundaries 
of vehicle attitude changes, but this procedure is strictly up to the 
discretion of the user. 

A note of mathematical caution about the use of correction vents: no single 
vent can account for the effect of two distinct vents on a vehicle 
trajectory. For example, consider the following simplified timeline. 



Suppose the vehicle experiences two distinct vents, v^ and v 2 , indicated 
above by the solid lines. Mathematically, there is no single average vent, 
i.e., dotted line, which will produce the same vehicle trajectory after 
time, T, as two distinct vents, v^ and v 2 . See Reference 3 for a further 
discussion of this idea. The attempt to attribute many unknown error 
vent sources to one correction vent is, at best, an approximation. 
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3.1 THE VARIATIONAL EQUATIONS - INTEGRATION AND BOUNDS 

Consider the following simplified vent timeline with one correction vent. 



Vi 



r 

. i 



■Tii 

to tj t; 

n 

i t3 U 


v 

Attitude 1 

Attitude 2 



Note that the above timeline contains a vent straddling an attitude change. 
This is done for illustrative purposes only and is not the form one would 
implement on the computer. 

The equations of motion can be written as follows: 


'* * k k 

X = ? (x, x, t) + — Ao$o(t) ei + — A 0 $i(t) (vx + ex) 

M M 


k k 

+ — Ax(t) S 2 (t) (ex + v 2 ) + — Ax(t) s 3 (t) ex- 
M M 


(1) 


The function, ? (x, x, t), represents the accelerations due to all sources 
other than venting. The matrices, Aq and Ax(t), are transformation matrices 
from the body axis frame to the Mean' of 1950 Coordinate System (M0F50). 

The step functions (sj(t), i = 0,...,3> are defined by 

si(t) =1 t € (t 1# t i+1 ) 

s^(t) * 0 elsewhere. 


M is the vehicle mass (slugs) and k is a scale factor which converts ft/sec* 
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to Erad/;irin2 (earth radii per minute squared). Equation (1) Is in a 
Cartesian Mean of 1950 coordinate system with vents appearing in the body 
axis frame 1r. units of pounds force (lbf). 


Let the components of ej = 



and take the partial of Equation (1) 


with respect to parameter, zj. Then 


3x 


Note the presence of step discontinuities in Equation (2), vhe variational 
equation. Care must be taken to integrate the variational equations accurately. 
The integration of the state equations also involves integrating through 
vent force discontinuities, but the vents are a small contribution to 
the total acceleration experienced by the vehicle. In Equation (2), 
however, the discontinuities are the sole driving functions. Section 
4.0 contains a discussion of the effect of integration accuracy on the 
propagation of covariance matrices. A smoothing' scheme, similar to the 
one used for the vent model in HOPE, is required to properly Integrate 
the variational equations. 


3 1 3x 3 ? 3x k 

- * + -y — + --- A c s 0 (t) 

1 3 X 3zi 3 X 3zi M 


l\ k 

0 1+ — AgsiCt) 
.0/ M 


+ — Aj(t) S2(t) ( 0 ) + — Ai(t) S3(t) (o 
M \0/ M \ 0, 


( 2 ) 


Figure 3-1 contains a flowchart and variable description of the smoothing 
logic Implemented in integrating the variational equations. 


An additional comment concerning the figure on page 6 is that the correction 
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FIGURE 3-1 (Concluded) 

SMOOTHING LOGIC VARIATIONAL EQUATIONS 
FLOWCHART VARIABLES 

PART(1 )-PART(3) - forcing function in the right hand side of the 
variational equations. 

TBL0CK(1) - current integration time points 

* 

TBL0CK(4) - integration step size 
TTEST - integration averaging interval 
MASC0N(IDUM4) - start time of correction vent 
MASC0N(IDUM5) - stop time of correction vent 
RBECI - transformation matrix from body to M0F5O 

7 

K q - scaling constant which takes lb force to Erad/min 
W - vehicle weight 
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vent, %\ % is shown straddling an attitude change. Because of the implementation 
of the attitude routine, BODATT, in HOPE, attitude maneuvers are assumed 
to be instantaneous. This introduces additional discontinuities at points 
of attitude transition. The user would be advised to create two distinct 
correction vents out of the single vent, ej. 

Now consider the following simplified vent timeline: 

ei 


to ti 

E* LVLH *3 


The equations of motion are: 


x * ? (x, v, t) + — A (x, v) so(t) Rq 

M 



where v * x, A (x, v) is the transformation matrix which goes from the 
Local Vertical/Local Horizontal (LVLH) attitude to the Mean of 1950 system, 
and Rq is the constant transformation matrix which goes from the body axis 
frame to LVLH. The matrix, A (x, v), can be written as: 

.+ *►« + ♦ + ♦ 

+ + (x X v) X X • -X X v • -X 

A (x, v) = "> 7+ *7 7+7 

|x x v| |x| • jx X v| • |x| 



Oefine the vectors: 


a l 


/+ +. + 
(x X VJ X X 


-xxv 


a 2 = 


jx x v| 


and 


-x 


a 3 = 


i*i 


Rewrite Equation (3) as: 


x = f (x, v, t) + — (c* ti (x, v) + c 2 a 2 (x, v) + c 3 a 3 (t, v)), 

M 



3v k 3q 3c 2 . 3c 3 . 

— + + — (-- q + --- a 2 + --- t 2 ) 

32^ 3x 3zj 3v 3zj M 3zj 3z^ 3z^ 


k 3ai 3x 3ai 3v 3ao 3x 3ao 3v 

+ — (q ( :J - — ♦ --- — ) ♦ c 2 — + --- —) 

3x 3zi 3v 3zi 


M 


3x 3zi 3y 3zj 


3a 3 3x 3 a -j 3v 

♦ C 3 (-*- + --f — )) . 

3x 3zi 3y 3zj 


( 5 ) 


The terms, 
and 


k 

(ci 

3q 


da 2 


3t 3 

3x 

M 


+ c 2 


c 3 

-r) 

— 


3x 


3x 

3x 

3 Z1 

k 

(ci 

3q 


3a 2 


3a 3 

3v 

— 


+ c 2 


c 3 


— 

M 


3y 


3v 

37 

3 Z1 


( 6 ) 


( 7 ) 
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3f 3x 3f 3v 

can be grouped, respectively, with -- and -- . 

3x 3zj 3v 3zj 


We shall now show that Terms (6) and (7) can be neglected in the variational 

+ ->• 

equations. It can be seen by examining the parti als of A (x, v) 

(Reference 4) that the following bounds exist: 


ii 


.+ 

3aj 

3x 


ii 







ii 






ii 



ii 




NOTE: || • || indicates matrix norm. 


(8) 


The constants, lq, k 2 , k 3 , k 4 , depend on the eccentricity of the orbit. 
For small eccentricities the constants are very close to unity. 


Consider the following table which compares the relative magnitude of 
accelerations affecting the Shuttle with a weight of 200,000 lbs. 


120 NM Circular Orbit 


Acceleration Source 
Earth central body 
Earth perturbations ( 02 ) 
Venting (40 Ibf) 

Drag 

Sun-Moon perturbations 

Erad = 


Acceleration 
18 Erad/hr 2 

- .02 Erad/hr 2 

4.0 x 10* 3 Erad/hr 2 

- 3.0 x 10" 4 Erad/hr 2 

- 1.0 x 10“ 6 Erad/hr 2 
Earth radii 


at ax a? av 

Let us now examine the tern, -- + , which appears in the right 

3x azj dv az^ 


hand side of Equation (5). We shall only consider the contribution 
of the central force term to the total acceleration, f. 

The partial of the central gravity term is as follows: 


9 yx ]i 

I - 3y ~s.“c * 

3x |x| 3 Jx| 3 |x| 5 


-*->T 

xx 


(9) 


where p * 19.909 Erad 3 /hr 2 . If jxj 35 1.03 Erad, then the partial of the 
central gravity term = 36.4/hr 2 . 


Expression (6) can be bounded by 





where use has been made of bounds in Equation (8). Expression (7) can 
be bounded by 


3aj 3a2 da 3 1 

C 1 + c 2 + c 3 T* i II c II max ( k 3» k 4 ) --- 
dv dv 3v |v| 


(ID 


using the bounds in Equation (8). 


For a 40 lbf vent and a 120 nm orbit. Equation (10) can be further 
bounded by 3.8 x lQ-fyhr^. Thus, the gravity term (Expression (9)) is 
roughly 10 4 times more significant than Term (6) in the variational 
equations. To further bound Expression (11), we must make use of Kepler's 
Third Law for the period of a circular orbit. For a circuHr orbit, 

i ; 



( 12 ) 


But we also know that T 


so that 




2» |x| 

"\W 


Hence, Expression (7) can be bounded by 


3a| da2 da3 

C 1 T + c 2 --- + c 3 --- < || c || max (k 3 , k 4 ) 
dv 3v dv 


1 ^ 11 C 11 |J|l/2 

„ < .... 


or 

< 1.4 x 10" 2 /hr 3 (13) 
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Gravity (Expression (9)) is roughly 103 times more significant than 
this term. 

Thus it appears that for orbits of small ellipticity, terms of the form 
aA ^ . * A ft. v) 

Jg" ' x » V J ““■*“+ , represent a relatively insignificant contribu- 

tion to the variational equations. These terms have not been included in the 

HOPE/VENT program, either in solving for venting parameters or solving for any 
other dynamic parameter. 
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3.2 PROGRAM MODIFICATION GUIDELINES 

The guiding philosophy in implementing solve/consider venting forces in the 
HOPE/VENT program has been to use as much of the original code as possible. 

A little used capability to solve or consider Lunar MASCON parameters has 
been modified to accept venting parameters. Thus, instead of creating 
separate routines in the input processor, trajectory, and DC links of the 
program, the MASCON capability was used to set appropriate flags, dimension 
variables, and create the necessary memory locations.' At the appropriate 
stage, however, MASCON variational equations were replaced with venting 
parameter variational equations. 

The reasons for selecting MASCON parameters are threefold. First, the 
venting parameters are dynamic variables which must be integrated along a 
trajectory. The MASCON parameters are dynamic and thus no changes in the 
trajectory link would have to be made to accommodate new dynamic solve for 
or consider variables. Second, the venting parameters are an expandable 
set of variables depending, in number, on the program user. The MASCONs 
are also an expandable set of solve for or consider variables with a maximum 
of 100 MASCONS available to be solved for or considered. Thus, the requirements 
for the venting force parameters are matched perfectly by the MASCON parameter 
capability. A final, compelling reason for selecting existing code is 
to limit the size of the HOPE program. The trajectory link is the critical 
program link In this respect, with future trajectory modifications of 
HOPE being severely restricted in Size. There are currently only 4,000 
words of unused core remaining in the TRAJ link with which to add subroutines 
or improve present capability. 


3.3 COOE MODIFICATIONS 

The following pages describe which subroutines have been modified in HOPE/MDC 
2.0 to create the HOPE/VENT program. Also described are the extent of 
the changes and the reasons for the modifications. For partial listings of the 
modified routines, see Tables B-l through B-9 in Appendix B. 

Lines of 


Subroutine 

Code 

Comments 

ASSIGN 

2 

Alter pointers in VSTR array for trajectory link to 
allow MASCON and Earth gravity parameters in core 
simultaneously. 

DAUX 

17 

Stores off initial input value of VTAtJ array and 
calls MASACC when correction vent capability 
requested. 

MASACC 

72 

•Computes acceleration of correction vents in the 
state and variational equations. 

MCNPRC 

2 

Writes MASCON information which has been processed by 
the input processor on drum. 

TRAJ 

14 

Resets the value of the VTAB array to its initial 
value upon exiting TRAJ routine. 

TRIGER 

4 

Removes code which would branch around MASCON 
capability if the central body is not the Moon. 

TRAJRD 

2 

Reads trajectory information from drum into VSTR 
array. Changes allow TRAJRD to read Earth gravity 
parameters and MASCON parameters simultaneously. 

TRJSUP 

2 

Supervisory routine which initializes VSTR array for 
routines in TRAJ link. Changed TRJSUP to initialize 
MASCON variables when Moon is not the central 
body. 

HOPE 

4 

Change the map to reduce program run time. 
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4.0 RESULTS 


Program verification was accomplished in two distinct phases. First the 
solve for capability was tested using the HOPE dummy data capability. The 
current venting logic in HOPE/MDC 2.0 was used to generate a trajectory data 
tape which contained known vents. HOPE/VENT was then used to solve for the 
known vents. Secondly, the consider capability was tested by comparing 
answers with the IGS Burn model in HOPE. The consider test phase was a more 
Involved one because the 16S Burn model and the venting logic use different 
methods to achieve the same goal. The differences in the way the 
trajectories are integrated is a prime example. 

In addition to the above verification efforts, other runs were made to test 
program flags and pointers. The HOPE program was run with a checkmode print 
of the input processor, and program flow and flags were examined. 

The pointers for the VSTR array were examined to ensure that using MASCON 
subroutines to solve for vent forces does not present hidden pitfalls. 

Special attention was given to routines which either read or write MASCON 
information on a drum. The MCNPRC and TRAJRD programs were examined to ensure 
proper input-output of MASCON parameters. Using COMGEN CCREF, all variables 
which relate to the proper functioning of the MASCON routines were traced to 
their occurrence in the HOPE code and studied to see if their occurrence 
would cause any problems in solving for or considering vent forces. The input, 
trajectory, differential correction (DC), and covariance propagation links of 
the program were dumped using checkmode print. In addition, snapshot core 
dumps of HOPE were made during various links to determine whether key flags 
had the correct values. 
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Two stand-alone checkcases were also employed for verification ("stand- alone” 
in the sense that no external comparisons were required). The two cases 
studied were: 

1. Doing a DC run starting with the exact values of the solve for variables. 

2. Closure tests on the integration of the vent force variational equations. 

Test #1 produced excellent trajectory agreement between the dummy data and DC 
runs (to fifteen digits). Exact agreement was not obtained because 
observation computations are supervised by different routines in the DUMDAT 
and DC processors. OBCOMP is the supervisor in the DC mode while DUMCAL 
supervises in the DUMDAT mode, and differences in the order of the 
arithmetical computations give slightly different answers. 

The second test was designed to examine how discontinuities in the varia- 
tional equations affected integration performance. At first, discontinuous 
forcing functions in the variational equations were evaluated at their exact 
values. This was found to be inadequate for the propagation of covariance 
matrices and a smoothing scheme had to be adopted in the variational 
equations, similar to the scheme for smoothing vent discontinuities in the 
state equations. Performance of the integrator improved considerably. For 
further discussion of this point, see the results of the consider test case 
(Section 4.2). See also Figure 3-1 for a flowchart of the variational 
equations smoothing logic. 

4.1 SOLVE FOR VARIABLES 

Dummy data with known vents were generated using HOPE/MDC 2.0, Two distinct 
tapes were created, one with noisy data and one with perfect data. The 
perfect data were used as a checkcase for one of the stand-alone tests of the 
HOPE/VENT program. The on-off times of the vents were input to r^.'F/VENT 
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and the program was executed to solve for the unknown vents as well as other 
selected parameters. 

In Case I, the program was executed to solve for both the state and vents. 

The state was perturbed by 500 meters and the program was required to 

solve for a 20 % vent. Perfect data were used in Case I. Due to the 

formulation of the HOPE iteration stopping rule, however, perfect data 

will cause more iterations than are normal. This explains why seven iterations 

were required in Case I. Case I results are given in Table 4-1. 

In Case II, only vents were solved for and again, perfect data were used. 

The final iteration produced deltas in the solved for variables of 10" 4 lb f , 
clearly an insignificant correction. Case II results are given in Table 4-2. 

In Case III (Table 4-3), vents and drag were solved for and the state was 

fixed at the exact, true value. The program was able to determine drag to 
within 98* of the true value. 

In Case IV (Table 4-4), two contiguous vents on the timeline were solved for 
and the program did an excellent job in the solution for each. 

Finally, in Case V (Table 4-5), SWCO data were used to test the program in a 
non-HOPE generated environment. The vent being solved for was a small vent, 
extending for 4-1/2 hrs. The program was not very successful at solving for 
this small vent, making an error of 50* in magnitude. Detailed studies will 
have to be performed to determine the feasibility of solving for small vents. 

In all the cases in this section, the vehicle attitude was Inertial Hold. 
Additional studies will have to be performed using different attitude modes. 
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TABLE 4-1.- CONVERGENCE RESULTS - SOLVE FOR CAPABILITY 


Case I - Solve for State and Vent 




Iteration 



True 


#0 

(Initial) 

#1 

#3 

#5 

#7 

(final) 

Value 

X 

-2806400. 

-2806399.8 

-2805968.7 

-2806073.3 

-2806083.4 

-2806084.5 

Y 

2878400. 

2876877.5 

2878177.4 

2878115.3 

2878110.3 

2878109.7 

Z 

5259600. 

5259190.3 

5259148.9 

5259323.7 

5259340.4 

5259342.2 

X 

• 

-7.04 

-7.7035 

-7.188 

-7.054 

-7.0413 

-7.04 

Y 

• 

-6810. 

-6810.19 

-6809.98 

-6809.99 

-6810.00 

-6810. 

Z 

3720. 

3720.13 

3720.32 

3720.03 

3720. 

3720. 

VENT1-X 

0.0 

-1.7369 

-.095 

-.0065 

-.00059 

0.0 

VENT2-Y 

0.0 

-8.9178 

-19.125 

-19.915 

-19.99 

-20.0 

VENT3-Z 

0.0 

.0267 

-.469 

-.0437 

-.00416 

0.0 

State Coordinate System: M0F50 Cartesian 
Vent Coordinate System: Body 

State - meters, seconds; Vents-lb force 
State RSS Position Error (Final - True): 1.99 1 

Vehicle 

meters 

Attitude - Inertial Hold 


TABLE 4-2.- CONVERGENCE RESULTS 

- SOLVE FOR CAPABILITY 

* 


Case II - Solve for Vent, State Fixed at True Value 

, 



iter 

'at ion 



Irue 


Ml 

#1 

#3 - 

. #5 

#7 

(final) 

Value 

VENT1-X 

0.0 

-1.717 

-.01758 

.804 0-5 

.113 0-8 

0.0 

VENT2-Y 

0.0 

-6.874 

-20.193 

■20.00 

-20.0 

-20. 

VENT3-Z 

0.0 

7.709 

-.0168 

.1196 0-4 

.2628 0-8 

0.0 


Units: Vents - lb force 


Vehicle Attitude - Inertial Hold 
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TABLE 4-3.- CONVERGENCE RESULTS - SOLVE FOR CAPABILITY 


Case III - Solve for Vent and Drag, State Fixed at True Value 





"Iteration 



frue 


10 

(initial) 

#1 

#2 

#3 

#5 

(final) 

Value 

VENT1-X 

0.0 

-2.384 

.01588 

.02249 

.23 

0.0 

VENT1-Y 

0.0 

-18.8 

-19.989 

-19.99 

-19.99 

-20.0 

VENT1-Z 

0.0 

.02976 

.611 D-3 

.012 

.0127 

0.0 

DRAG 

500 

3889 

1385 

995 

• 

983 

1000. 

Units: 

Vents - lb force 
Drag - Area (ft 2 ) 


Vehicle Attitude - 

Inertial Hold 


. TABLE 4-4.- CONVERGENCE RESULTS 

- SOLVE FOR CAPABILITY 



Case IV - Solve for 2 Vents, State Fixed 

at True Value 





iteration 



True ~ 


#0 

(initial) 

#1 

#2 

#3 

#4 

(final) 

Value 

VENT1-X 

0.0 

-3.23 

.065 

.0285 

.0285 

0.0 

VENT1-Y 

0.0 

-46.54 

' 50.035 

-49.99 

-49.99 

-50.0 

VENT1-Z 

0.0 

1.66 

.0607 

.011 

.011 

0.0 

VENT2-X 

0.0 

-39.89 

-30.01 

-29.97 

-26.975 

-30. 

VENT2-Y 

0.0 

-18.75 

-20.164 

-19.99 

-19.996 

-20. 

VENT2-Z 

0.0 

-6.705 

-.0417 

.0269 

.0269 

0.0 


Units: Vents - lb force Vehicle Attitude - Inertial Hold 
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TABLE 4-5.- CONVERGENCE RESULTS - SOLVE FOR CAPABILITY 
Case V - Solve for Vent, SWCO Data Tape, State Fixed at True Value 

" " " ”” " “ : True 


fteratTon 



to 

(initial) 

#1 

#2 

(final) 

vaiuc 

VENT1-X 

0.0 

-.1335 

-.13586 

0.0 

VENT1-Y 

0.0 

-.4936 

-.4939 

-.5 

VENT1-Z 

0.0 

.1805 

.1305 

0.0 

Units: 

Vents - lb force 

Vehicle Attitude - 

Inertial Hold 


4.2 CONSIDER VARIABLES 

The consider capability was tested by matching answers obtained with the 
H0PE/B8.02 program running in the IGS Burn mode. HOPE/VENT was first modified 
so that the DOPPLER, JACHIA and TRAJ subroutines agreed with the subroutines 
in B8.02. Initially, discontinuous functions in the variational equations 
were evaluated at their true values. However, because the HOPE integration 
accuracy criterion applies only to the state equations, the integrator 
in the auto mode could not respond adequately to jumps in the variational 
functions. The propagated covariance matrices were not in agreement between 
this initial form of the HOPE/VENT program and B8.02. To substantiate the 
idea that integrator inaccuracy was the cause of the lack of agreement, the 
apoch was moved to the center of the integration time period and closure 
results were observed. Closure of the variational equations was poor and 
errors in propagated covariance matrices could be attributed to the errors in 
the variational equations. The integration model was changed in order to 
smooth the variational discontinuities and immediate Improvement was observed 
In the results. 
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The following is a description of the consider test case: 

1. 5 min APU vent 

2. Inertial Hold Attitude 

3. Observations precede beginning of vent, i.e., 


To 
(epoch) 


mil 1 1 mi n mini mi ts s 

5 min APU 
Vent 


Covariance propagation times 

J LJ 1__L_ 


Observations 
300 Observation Frames 


4. Consider the uncertainties in APU vent components upon the state 
covariance matrix. 

5. Start with the exact state and do a one- iteration DC run. 

APU VENT DEFINITION 


Force, Body Coordinates 

Value 

Uncertainties 

F* 

2.85 lb f 

oF x - 1.43 lbf 

F y 

2.28 lb f 

<jFy a .52 lbf 

F Z 

17.17 lb f 

cF z s 7.25 lbf 


Tables 4-6 and 4-7 give the IGS Burn values of the state and covariance 
matrix at two comparison time points. Tables 4-8 and 4-9 give the HOPE/VENT 
values of the state and covariance matrix, with variational equation smoothing, 
at the same time points. At 2 hrs and 42 min after vent cutoff (April 1, 

20 hr, 0 min), one observes a 1.3 meter RSS difference between the IGS 
Burn and HOPE/VENT values of the state, and a maximum deviation in the 
position sigmas of .2 meter. At 6 hrs and 42 min after vent cutoff 
i-pril 2, 0 hr, 0 min), one observes a 3.6 meter RSS difference In state 
values and a maximum deviation in position sigmas of .5 meter. Of course. 
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since the variance is the square of the sigmas, the observed differences 
in the ^covari ance matrix along the diagonals is much greater. 

To account for the deviation in the propagated covariance matrices between 
I6S Burn and HOPE/VENT, the integrator was modified to function like the 
integrator ih the IGS Burn: model. The integrator will integrate to vent 
discontinuities, stop, and restart in the Runge-Kutta mode. This ensures 
that no difference table entry will ever straddle a vent boundary. The 
modified integrator was checked out on a rigorous vent timeline and the 
value of the vent force assigned at each time point was carefully examined. 
Tables 4-10 and 4-11 give the values of the state and covariance matrix 
for the modified integrator. The agreement with the IGS Burn model is 
excellent. The state position differences are on the order of 10~ 4 meters 
and the covariance matrices agree to all digits. 

Though the modified integrator is an unverified program change, it shows 
such marked superiority over the smoothing logic that the author feels 
it should eventually be adopted as the standard for venting. For the 
present, however, even though the smoothing logic did exhibit errors in 
propagating covariance matrices, the size of the errors does not preclude 
one from using this model. The author believes the errors in the smoothing 
scheme are acceptable and the first version of HOPE/VENT should be coded 
with smoothing logic for the vent force variational equations. 
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COVARIANCE PROPAGATION 
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5.0 USER'S GUIDE 


The following instructions refer to HOPE version MDC 3.0 which should 
be the first released version incorporating the capability to solve or 
consider vent forces. An update to the following instructions will be 
published in a working paper to accommodate cosmetic changes which are 
likely to occur in the variable names. 

As currently envisioned, the execution of the HOPE vent capability will 
require a separate map, distinct from the current HOPE map. This has 
been necessitated because of the desire to preserve the MASCON capability 
in its original form. The recently modified subroutines will be copied 
onto the PCF tape with changed element names. The original MASCON subroutines 
and the modified MASCON subroutines will exist on the PCF tape simultaneously 
and the HOPE map will determine which copy of the MASCON subroutines will 
be brought into core. 

The capability to solve for or consider vents must be used in conjunction 
with the current HOPE vent capability. The rule which must be strictly 
observed is that any end point of a correction vent must correspond to 
an end point of a nominal vent in the VTAB table. The simplest way to 
ensure that this will occur is to input, in the $VITVTL section, a nominal' 
vent with zeroes in all components which have the same start and stop 
times as the correction vent. If the user omits to enter any nominal 
vents while attempting to solve or consider correction vents, the program 
will terminate. 

The following is a list of additional strict guidelines that must be followed: 
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1. The program Is written for a maximum of 50 nominal entries in the 
VIAB table. If the VTAB array must exceed this size, the temporary 
storage array VOTAB will have to be re-dimensioned in OAUX and TRAJ. 

2. If HOPE/VENT is used* the scale card for the MASCONS must be input 

. as: SCLMC * ERAOI I, RADI AN, MINUTE* - 

3. The correction vent information is input in the $MASC0N section using 
the MC card, Ue., MCi * xj, x 2 , x 3 , X 4 , x 5 , x 6 , X 7 , x 8 , x 9 * 

The vari abl e, 1 , is the nurnuar of the correction vent; k \» x 2 » X 3 

• are, ^spectiyely, the body coordinate XYZ components of the correction 
vent in lbf; are zeroes; xg, X 7 are the start and stop times 

of the correction vent in internal time, i.e., minutes past midnight 
day of epoch; x 8 , xg are zeroes. The fields, xg and x 9 , may be omitted 
if desired. Fields *4 and xg must be input as zeroes. The correction 
vents should be listed in numerical order beginning with i = 1 and 
proceeding to i = N, where N is the maximum number of vents to be 
solved for or considered 

4. To solve for or consider a correction vent, the user must specify 
the variables, MASi, MLTi, MLGi, on the SOLVE or CONSID cards. The 
variable, MASi, indicates the x component of the 1th correction vent, 

MLTi indicates the y component of the ith vent, and MLGi indicates 
the z component of the ith vent. If more than one correction vent 

is being solved for or considered, the i-indices should be arranged 
In ascending numerical order on the SOLVE or CONSID cards. If correction 
vents are to be considered, the HOPE nominal ordering is used, with 
MASi, MLTi, MLGi representing the XYZ componerts of the ith correction 
vent in the consider covariance matrix. 
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The labels, SMASCON, MC, MAS MIT, and MLG, are temporary and are likely 
to be changed in HOPE version MOC 3.0. As mentiorsd, a working paper 
detailing the changes will be forthcoming. 

Since the current nominal venting subroutines are designed for use 
with only one vehicle, the correction vent logic has the same limitation. 
The correction vents must apply to the nominal vent timeline of vehicle 
number one. 


6.0 CONCLUSIONS AND RECOMMENDATIONS 

HOPE/VENT is a working modification to HOPE/MDC 2.0 which can solve or 
consider vent forces. The mathematical logic of HOPE/VENT has been verified 
and determined to function correctly. There exists only one reservation 
concerning the validity of the program and that is the integration of • 
the vent force variational equations. It is the recontnendation of this 
author that the cumbersome smoothing schemes in both the state and variational 
equations be replaced by an integrator which stops at discontinuities- 

i 

and restarts. Such an integrator already exists and has been partially 
tested by the author, but further verification would be necessary before 
it could become a part of a baseline HOPE version. 

Additional studies are planned, however, to fully test the performance 
of the HOPE/VENT program in its present configuration. The limitations 
and full capabilities of solving or considering vent forces will be explored 
in depth. 
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APPENDIX A 


MATHEMATICAL BACKGROUND 




which usually depicjts the state variables. Observations will be denoted as 
y 1 « F i (x, z) + n i, 

where the solve for parameters are indicated explicitly. In reality, 
however, the observations depend on the position and velocity of the vehicle 
at some time and the position and velocity may, themselves, be functions 
of the solve/consider variables. The term, r.j, is a noise term. The 
observation vector will be denoted as 

F 1 (x, z)\ in i 

r ; 

k F n (*» z) / n n 

where arrows used to denote vectors have been dropped for convenience and it 
Is understood that components of a vector have indices. Continuing with 
this representation, an actual observation will be denoted as 
y* * F (x, z) + n*. 
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Let xq, zq denote initial guesses for the solve/consider parameters. 
Define 2 ■ E ((xq - x) (xq - x) T ) . 

«r 

It can easily be shown that if y* is an actual observation and xq is an 
estimate for x, then the sequence of iterative values which minimizes 
the weighted squares of residuals is 


x 1+1 = Xi + (H + Ai T W A,)” 1 (Ai T W (y* - F (xj, z 0 ) 
+ 2" 1 (x 0 - X*)), 

where 


I ••• * • • • 


3F * 

Ai * (*i, Z 0 ) = 

dX 


3xx 


3F, 


<* x n \ 


(A.l) 


3F. 


Van •'•"ax-J 


(X-j* z q) 


Let x* denote the solution to the minimization problem. Then x* satisfies 
the equation 

A* t W (y* - F (x*. zq)) + 2' 1 (xq - x*) * 0. (A.2) 


Suppose that z 0 is perturbed slightly to the new value, z 0 + Az. How 
will the optimal solution, x*, change? If Az is sufficiently small, it 
can be shown that 

Ax* » (r 1 + A* t W A*)’ 1 (A* t W Ay + 2 _1 (xq - x*)) , (A.3) 


where 

Ay* * y* - p (x*. z 0 + Az) 

and Ax* is the perturbation to the optimal solution, x*. Now 


A-3 


3F 

F (x*, Zq + Az) a F (x*, z 0 ) + — (x*, Zq) (A.4) 

3z 

so that 

dF 

ta* ••(!-! + A, t HA*)' 1 (A* T My* - F (x», zq)) - A, t M — (x*. z 0 ) Az 

dZ 

♦ rl (xo - x*)) = - **A* t M 8* Az , (A.S) 

where 

*, ■ (A, T u a*)' 1 
and 

3F 

6* = — (x*. Zq) . 

3z 

Equation (A. 2) was used to simplify equation (A. 5) in the next to last 
equality. In terms of partial derivatives. Equation (A.5) can be rewritten 
as 

* 3x* 

** A* t W B* . (A. 6)* 

dz 

Equation (A. 6) Is a sensitivity equation, telling which perturbation of 
a consider parameter would change the optimal solution the most. 

A more meaningful estimate of uncertainties in the solve for parameters 

due to uncertainties in the consider parameters is provided by the covariance 

matrix. 

Writing, once again, the iterative scheme for the optimal solution: 


P 
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X i+1 * X) + (z- 1 + A, 1 H A))- 1 lA, T M (y* - F (x,, zq)) 
. ♦ Z' 1 (xj - x<» . 


(A.7) 


Let x and z be the true values of the variables being estimated. Define 

A 

ei+1 * x 1+1 - x . 

Now 

Ay i * y* - F (x 1f z 0 ) = n + F (x, z) - F (x is z 0 ) * 
n - Aj + Bj Az, where Az * z - zq . 

Subtracting x from both sides of Equation (A.7) gives 

"•1 A / A 

e 1+l = e i + ’h (-A^ W A-j e-j + A-j^ W Bj Az + A^ W n + £ (x 0 - x^). 


where 

= Ui T M A, + 2” 1 )" 1 . 


(A. 8) 


Then 

e i+l * ^1 (“A-j"** 14 Ai + A^ W Bj Az + A^ 14 n 

+ E“*> (xq - xj) + Aj T U Ai (A. 9) 

» (Ai t W Bj Az + Ai T U n + E’ 1 (So - x)J 


E (ei+l) * 4»i (A^ WBi E (Az) + Af T 14 E (n) + 2“* E (xq - x)) 

Assuming that 

a) E (Az) * 0 , 

b) E (n) ■ 0 , 

c) E (xq - x) ■ 0, 


then 

E( c^ + ]) * 0, and we have an unbiased estimator of x. 
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Forming the outer product* of Equation (A.9) with itself and taking the 
expectation of both sides gives 


E C e t+l c i^+l) * .'♦i ^ A i T H B i E Az B i T w A i 
; + A t + 2‘V + Aj^ W Bi Eaz»Ax 0 2 * 

+ ^Axo.Az B i T « Ai) , 


where 


(A. 10) 


*to 0 ,ta T * I te,te ; 0 '" E * *0> «So'- *) T J 

E&z = E C(z - zg) {z - z 0 ) T ) « 

Recognizing that 

* 1 * (Ai T W A, + r 1 )- 1 , 

we can rewrite Equation (A. 10) evaluated at the converged solution, x*. as 

C 0 • + ♦* U* T U B* Bj K A, * aJ W B, E 4z ,Axo s '* (ft u) 

+ £”^ 2 ^x 0 , Az ^ A *^** 

Here we have made use of the definition, Cg * E (e-j+1 e i+l^)* 

Irfe can also obtain the cross covariance matrix, E (e i+ l (z - zq) T ). 

If we multiply equation (A.9) by z - zg and take the expectation of 

both sides, we obtain 

CRg - E Ui+1 ( z - *0) T ) a ♦* ( a * T W b * E Az + E " 1 e Ax 0 ,Az) ( A ‘ 12 ) 

* The outer product of a column vector, x, is xx^. 
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Suppose we are solving for the state at time, t 0 . Then Equations (A. 11) 
and (A. 12) contain submatrices which contain the covariance of the state 
and the cross covariance of the state with all other solve/ consider parameters. 


In order to derive the equations for the propagation of the state covariance, 
we must depart from previous notation. Let x(t) * Ui(t), X2(t),...,xg(t)3 
> denote the state and let p = (pi,...,p m ) denote the set of all solve and 
consider parameters minus the state variables. Let the state propagation 
equation be written in simplified form as 

x(t) = f (x (t 0 ), p, t) (A. 13) 


A A 

If x (t 0 ), p are our best estimates of the true values, x(t 0 ), p, then 
x (t) = f (x (t 0 ), p, t) , . (A. 14) 


and 

A A 3X(t) A 

f (X (t 0 ), p, t) = f (x(t 0 ), p, t) + (x.(t 0 ) - X (t 0 )) 

3x(t 0 ) 

3x(t) 

+ (p - p) . 

3p 


We can rewrite equation (A. 15) as 


a 3x(t) a 3x(t) 

x(t) - x(t) (x(t 0 ) - x(t 0 ) + (p - p) . 

3x(t 0 ) dp 

(A. 16) 

Denote A^t) = E [(x(t) - x(t)) (x(t) - x(t)) T ) , 

(A. 17) 

ip 

A Ax,Ap (*) a E CCx(t) - x(t)) (p - p) ), 

(A. 18) 

and a a _ 

A Ax,Ap * E C(p - p) (p - P) T 3 . 

(A. 19) 


Taking the outer product of Equation (A. 16) with itself we obtain 


A-7 


(A.20X 


A Ax(t) 


3x(t) 


3x(t)T 3 x (t) 


’»*'(‘o) A4x (to) Tt'd'S IxnjV^P (t o) • 


5x(t) ' 3x(t) _ 3x(t) T 3x(t) 3 x (t)T 


— + 

3p dp 


A tvx,Ap T (to) + A Ap,Ap 


3x(t d ) 3p 


3p 


The cross covariance becomes 


3x(t) 


3x(t) 


V" (t) . ------ AlMp (t o) ft4p>4p . 


(A. 21} 


This completes the mathematical background for sol 


ve or consider parameters, 


APPEND I a 8 

PARTIAL SUBROUTINE LISTINGS 



B-I 


The following listings of modified subroutines list only the general 
area of code affected and are not complete representations of the altered 
subroutines. For a complete listing of the modified code as well 
as the original code (unmodified H0PE/M9C 2.0), please contact the 
author. 


Table Title Page 

B-l Partial Listing Subroutine ASSIGN B-l.l 

B-2 Partial Listing Subroutine DAUX B-2,1 

B-3 Partial Listing Subroutine MASACC B-3,1 

B-4 Partial Listing Subroutine MCNPRC B-4,1 

B-5 Partial Listing Subroutine TRAJ.. B-5J 

B-6 Partial Listing Subroutine TRIGER B-6.1 

B-7 Partial Listing Subroutine TRAJRD B-7.1 

B-8 Partial Listing Subroutine TRJSUP B-8,1 

B-9 Partial Listing Subroutine HOPE B-9.1 


/ 


Table B-l - Partial Listing Subroutine ASSIGN 


CD 

I 


9MS? 

S9153 

e?isi 

2015 S 
9015#. 
22157 
221*0 
9141 

r?i 6 ? 

221*3 

9«145 

9016* 

90**7 

9?i70 

9?17 J 

*1177 

90173 

9C171 

99175 

9?17* 

99^7 * 

99177 

C?500 

2?5r,o 

92517 

9251 3 

295| A 

29520 

99573 

92523 

29573 

29523- 

225 25 

2952* 

22527 

'*2530 

*'953! 

29537 

22533 

92531 

92535 

C953A 

*9537 

92510 

2251 1 

92517 

29513 

99511 

2051* 

92517 

22550 
22557 
22553 

22551 
29555 
2255 A 
295*0 


* 


*;*• 

*♦7* 

**«• 

*99* 

700* 

70)6 

707* 

7036 

7 Pi* 
70? • 
70*6 
7076 
70 A • 
709* 
71*. 
7t I* 
71?* 

71.. 
7116 

7 I 5* 
7U* 
717# 
7106 
7196 
7706 
7716 
7 7?# 
723- 
7?1# 
725 . 
77&. 
7276 
7 2,1 . 
7296 
7 3r>. 

7 3 | . 

7 3 7. 

73.. 

7 3 1 • 
7356 
736. 
737# 
73a. 
739. 
71a. 
71 1 . 
7176 
713. 
7li« 
715# 
71*. 
717 . 
711. 
719. 
75n. 
75 1 • 
757. 
7536 


Kl 1*1 
M 13) 
MID 
K( 171 
K( I A) 
" K * 19) 
K 1 2 f* ) 
* ( 21 I 
M2? ) 
Rf 231 
K121I 
K 1 25 1 
K I 26 1 
K < 27) 
K ( ?B 1 
M 29 ) 
K (30) 
K t 33) 
K ( 3<j » 
K< 31 ) 
K(32) 


I 


"AUG** . 

hath* l • 

«STA*I 

MGOS! - * 

ML A MO 1 

MG1M . 

MG A A 1 
MOAA I 
moha 1 
MST*«I 
MGB 1 ) 

M0B5T 
MOA 5 I 
MOfi| P T 

moclpt 

MMP4 VC 
M A y v A 1 

MMoA VC* NS ° LVf * 1 * ,50 ‘- VE * 
MAVECl . | 


I 


• 1 
• 1 
1 I 
• I 
- ! 


I > • I 


• Noei>(i«) i 

W.7ITE ( KDRUM2 ) j {PojNTfj I » | • 

’cALi ASSCKM 12 , 2* ) * * * * 1 " 

lio 00 150 1 • |7 , *o * 

159 IPOiNT(f) . 0 

00 1*6 I m I , 1 K 
1*9 Kll) • O ’ U 

| ♦5l!6N F tT-APpi. V-1T5UM VAP1A0LC STORAGE 


• LpOJNT 

♦ LK 


: v 


17' 


MAMpLV ■'MSTARI 
DATA ! I * MAPOLT 



NSOLVE ♦ MVC(I'I) ♦ NVC1I.2) 
NVCi?,!) ♦ NVC<2.2> * 

* “ ‘ ) 

) 


( 2 , 1 ) 

( NAtRAT 
I NSOLVE 
NSOLVP 
( NSPLVE 
( N501.VE 


N50LVE 

NCMST5 


NCMSTS 

MCNSTS 




NCM5TR • NCMSTS ♦ 

. MOLOO ■ MatATI ♦ 7 • 

MOP » HOLO 6 ♦ 2 • 

MLABL5 “ MOO 626 

MSCALS » MLARLS ♦ 2 • 

HANOI ■ Msf A l5 6 2* , .... 

MC p R3 » MONOS ♦ 2*NS0LVF 

H1GH3 • MUPB3 * .* 6 ° U NLP11, 

ML**R1 ■ MIr,Al • SfJ*Mf6(ll 
M1GH1 . M1.PH1 * 1* r'iipu, 

MC* £J.1«M I Grm+59«N (G ( 2 ) . 

MCAFC3 « MCivEja 

IJlfMP ■ NJc?/0(l, + NCSrRO(i) 

* » jc *» ■«*,» ♦ > 

ITCHP. M A X 0 ( N J ( • ) t I I 

M C«Pl N ^ < H > i r t NC, !ri 'I** 0 1 ,T t MP ■ 0 

. ncSEi » M.JP3 ♦ »VEMp*> 

MCPHv»30MCSf3*NCS( I ).2 
DC6MC3 * Mf«MJ3 
f TEMP • NJcnO(?i ♦ NrirwO (7) 

*HE* ©I MCWMC3 » MCPMC3 ♦ 

HJH3 - HCRMC3 ♦ NCSCWf>(2> 


NUCPOif) ♦ 3 


• NE- 
66-1 


• NED 


...W 


'*'■‘5*1 


7 *ii 


JTfMP« MAXDNJf •*) , ( > 




Table B-l - Partial Listing Subroutine ASSIGN (Continued) 


C9 

« 

ro 


C35AI 

cosa* 

C05AH 

--•gag 

3.-5AA 


:•?! 67 

33**7 


33573 
33573 
3357*1 
33575 
33*0? 
PSAO? 
33 AO ? 
3" Ar; 7 
33*0 » 


""toi 

30505 


C360* 

33507 

3«A | 0 

33 A | | 

335 I 7 

rc a » 1 

3" 6 I *1 

335*5 

?3AI* 

33517 

3"A20 

3352* 

C 3 A 7 7 

33573 

3362*1 

33625 

33 *2A' 

C?677 

33630 

33631 

33632 
3" A 3 3 
3 3 A 3 *• 
3" A 35 
33636 
33537 
33 * HO 
C361 1 
336*12 
33**3 

C" A *M 


C3*HH 
3?6H5 
306 HA 
306HA 
3?6*0 
33*61 
C366H 
33**6 
03*71 


75«* 
755* 
75 A* 
757* 
7*n* 
753* 
7fn* 
7*|* 
76?« 
76^* 
7* H * 
7A5* 
7o A* 
767* 
7 A a • 
7 *?• 


I S 




ITEHP- MAX«fNJt?l.l ) 

IF < N J ( 2 * * NCS* 2 * •€ o. 0 » !TC«P ■ 0 
HCSH3 • Hjm.3 ♦ T t E h P** 

MMf 0N3 *mCSM 9 + DC 5 < 2 ) •? 

HMCC«3 ■ Mrif0Jl3 ♦ 2 • NPARR • NHCON _ , 

M£HOfA ■ KAX-fM f SE3 ♦ 2*NCS(|) t MCSH3 ♦ 2*NCS»2>» 
• MHfr.W3iH!0BH "* S?*N|<;<2> ) • t 
• S30.500»Hfl0 


IP (MEnDFA-lVSTr) 
Hr** IFATAL ■ IfaTAL ♦ 
CALL FLIP* * ,!FLpJ 
*H1TC iF0OT,70Ce* 
60 TO 5H0 


IlINKlS)* MENOf*. LVSTP 


772* 
773* 
77h* 
775* 
77a* 
777* 
7 7a* 
77?* 
78f»* 

78 I • 
787* 
753* 
7R><* 
7H5* 
7 H A • 
7«7* 
78«* 
789* 
79n* 

79 J • 
79?. 
791* 
7 9 H • 
79$* 
79a* 

797* 

798* 

79?* 

809* 

801* 

8C7* 

903* 

93** 

805* 

8DA* 

007* 

809* 

809* 

Bln* 

8I|* 


COMPUTE REAO trjOICpS FoR F I T-arPLY-I TSUM 

I 

5 ** K i 1 ) » M OtO^ 

■ “ 1 
MUfn I 0 . 


fl 




X < 91 
Kl H ) 
K I 51 

XI 61 
XI 71 
K ( 6 » 
K f 9) 
X C 131 
X (11 J 
Kit?) 
K( I 3 » 
KMH) 

X I t 6 ) 
K ( t A) 

XI 17) 
Kllll 
X 1191 
KI7*)) 

X 12 I ) 

XI 2?) 
XI 231 
X I 2H 5 
X I 25 ) 
X I 7 6 ) 
X ( 7 7 ) 

X ( 28 ) 

XI ?9» 

X I 30 ) 

XI 31 ) 
X I 32 1 
X I 33 ) 
X I 3*4 ) 


MT»t 
Nt.P*3 
Hl6a3 
M I f, 8 3 
Ml PrtH 
Hi.pnH 
MJC,rH 
M t r, n H 
MlCiBH 
MC6F J3 
MCSF3 
MC«mJ3 
MCSm3 
MMeoM3 
MlCC*3 
MAppLY 

n* t a 1 1 

MAPP t 0 
MCNAI.3 
MAI.RN3 
MCHaLH-i 
Mai.eNH 
MATWAt - 
M L A M 0 | 

M G 8 H I - 
M G n 8 I 
M088I - 
M088 I 

hstaRI . I 

MTR1 

map? to - I 


- } 

» 1 

- 1 

♦ 50*M 1 6 ( 2 1 - I 

* 2*NCS<t> - * 

♦ 2*nCSI2> - I 

. 1 


- 1 


I 


I 


NDRUHIH91 I 


CALL RANLAS I KORUN? , . 

WRITE I X0RyM2 I j | Pof NT j J > * J - * 

*CALl AS5CXM I 3 1 '! I * 

5*» 00 650 I * | 1 I *01 NT 
SS M tPOINTItl *0 

00 5 AO T » | « lk 

Kill • 0 • 


. tpiM ; . 


•NEW 
• *-| 




y. 


* - 

1 


P?*73 

33673 

«* n l 1 t 


812* 

813* 


60 TO 709 

c assign vartadlf storage for dummy data 


V 










■MA***ai***MM*i*i 




L_ 


k 




Table B-l - Partial Listing Subroutine ASSIGN (Continued) 


10*9* 

ion* 

I01&* 

10*7* 

I P9A* 
1089* 
t 05*** 
1CS|* 
ICS?* 
loss* 
I05*»* 
less* 
j cs a* 

1C57* 

ms** 

ICS9* 
1C*"* 
I06|* 
106?* 
I 06 j • 
I C 6 * • 
I 065* 
IP66* 
1067* 

IPfcn* 
I C 6 9 • 
I 0 7*** 
1071* 
IC72* 
1073* 
107m* 
1 P7s* 
1076* 
1077* 
IC7»* 
1079* 
I OAn* 
»00|* 
1082* 
»0«d* 
IOPm* 

ioe=.* 

IC««* 

1087* 

teas* 

10B7* 
I0VP* 
I 09 | * 
107?# 
1093* 
I C?*»* 
ion* 

1096* 
1097* 
109** 
109** 
I 100* 
I 101* 


MSCC V* ■ MtBCVS ♦ 2*NPR 

HENOCR ■ M5CC9*» ♦ 2*NPR 
IF <MEnOCR-l VST p) 976 »976 »97S 
975 1FATAL • IFATAL ♦ I 

CALL FLlP<l,!FLp) , 

. WRITE (K0UT,2CC0» ILIN* «§» , MENOCR » LVSTR 
60 TO 981 


976 K{ »» 
KI 2) 
KI 3 I 

ki ** > 

• KI 5) 
K I A ) 
K l 7) 
KI 81 
KI 91 
K I t C I 


MAp. VI 
MAT* A* 
MAT** 4 * 
MAU6N7 
MAtr.f*7 
M A L G N 8 
MAlgH* 
MLBCV8 

ml nr V* 

MEHOCR 


CALL RANLAS<*0R||M2iN0RUM<591 » 

WRITE <K0RUM?> I I PO|nT < f • * I ■■ l 1 1.1*0 1 NT I 9 
* I K ( J 1 * I ■ l » LK > 

C AL| ASSCKM»8»I0» 

9*| 00 *02 I - I t LPOlNT 
982 IPOINTf I * a O 
00 988 !«1 .lK 
989 XU) • 0 

£ assign vehicle one trajectory VARIABLE storage 

986 NTM**2 - 2 • NYmAXI t) 

MYTR& ■ l 

€ MO I F I a My 786 ♦ 10 • NyTPG 

C • 

MYl * MfliFl ♦ 26 • NYMaXM) 

MYHI ■ MYl ♦ NYmA)<2 

MYPPI ■ MYP I ♦ NYMAX 2 

C MTTR6I ■ MYPP1 ♦ NYMAX2 . - 

C MAUXI » MTTR6I ♦ *» • NTTR(j|l» 

C MCNALI • M A (I X I ♦ 2 • NAUXt . -• 

MALCM • MCNALI ♦ 12 • MALI 

MCNl.Pl • MAlGNI ♦ 2« • NAL! 

MLNYI a MCMLPI ♦ 8 • NLP<ll 

MCNIGI a MLHTI ♦ MLPIlI -- - 

* ' M I N V | a HCmTGI ♦ 21 ♦ N f 6 I I I 

MLPBI • H I H Y t ♦ NlgU) 

M I 6R 1 a MLP81 * 16 * NLP C I I 

C ' MCPfJI « Mfr,81 * S3*NIG(1» 

MC»CCI a MCaEJI 

ITnMP - NJCWOUI ♦ NCSCWOII) ^ , 

IF 1 I TE M P *Nf* 3» McwEcl • MCWEC* ♦ N *>CRO||) *3 
MJFI a MCWECI ♦ NCsCwOd) 

I TEMPa MAX«<NJ( | > , I ) 

if i nj»»> ♦ Ncsni ,eq« o i Itenp » o 

HCSF I a Mjr| * jTEMp*2 
■” ' MCWMJ|a»iCSE|*2*NCSd 1 


$ 

%% 


i ir** 

I 1 P 3* 


MCWMC1 » NCtyMJI 

1 TEMP ■ NJCW012) ♦ NCSr*0<?> ^ - 

_»*. ... « . * »■• > . p • . Ma<.tKR I m UrWMP I * M.irwn # ? 1 ♦ 3 


k 


Table B-l - Partial Listing- Subroutine. ASSIGN (Continued) 


?|15I 

;iis3 

eus* 

01*5* 
Ol«S7 
51**0 
CM*? 
CIHA 3 
C|*4* 
S i 

civ* 

7 

51*70 

51*73 

' , l«i7<- 

5 1 *75 

51 * 7 * 
nsm 
ciso? 
5 1503 
515 IP 

2 is 1 1 

51517 

c 1517 

01517 
OlM* 
ms t<r 


1107* 
I IC)« 
110«l» 

1105* 
UP** 
1107* 
1 1 0 • • 
I lev* 

1 lift* 

tin* 

Ml 7 ft 

111 3ft 

1 I I ft ft 

1 I 1 5ft 

I I l Aft 
1117ft 

in?: 

U?5* 

nil* 

I I ?7ft 
1 1 23* 
1 1 7 ft • 

1125ft 


HC*HC 1 • HCftHJI 

HEMP • NJCWOI?) ft NCSrwOlft} 

*W ,T T«??«ci p l SfSSSJ,;,"'*"' 1 ♦ "•»'«!*> » > 



01515 

CI51A 

>1517 

-j ; 

cis?n 

‘ CD 

ClS? i 

t' if 

7152? 

CIS? 3 

r-W 

Cl 52* 

A. 

2|5 ?5 
c is?* 


Cl 527 


C 1 5,35 
C 1531 
51531 
01532 
51533 


5153ft 
51535 
5 ! 5 3 * 
51537 
5 1 590 
0 1 5ft 1 
?|5ft? 
515*3 
? 1 5 ft ft 
5| 5*5 
515*5 
ClS*7 
OJSSO 


1127* 

117ft* 

I 1 2 V* 

I ISPft v 
I I 3 1 * . ' 
Ml?* 
1131ft 
lilftft 

1135ft 
1134ft 
ii37ft ■; 
M3«* 

3 137ft 
1 t 4 - 7 * 
HU I* 

I 15 7 * 
115 3# 

I I ftftft 

1155* 
M1»* 
11*7* - 

11 ft 5 • 

I I *9* 

1 IE"# 

1 1 5 | • 

1 19 7 * 
llSlft 
US ft • 

1 155ft 
MS(« 
1157ft 
I ISPft 
1159ft 


I TEMP* MAX" ( VJt 9 | , | ) 

NJI 2 ) ft NCSI 21 .Fo. o I JTrHp ■ 0 

MCSMI ft M Jm ( ft fTCPpft? *• 

MMC0N1 ■ HC5**f ’♦ 2 ft NCS I 2 ) 

MMCCWIftMMCON I ♦ I 7 *NPAR 5 ) *NMC 0 M+ 7 #NMCON 

MVT ARftMMCCwiftNMrODE 

MRTAR ■ NtfTAft* *P»NVTft« 

MPT*P ft MBTftB ♦ 7 * • MB T A B 
- - " M ETAB * 8 *Nff AR 

KMC AG A* MjktaB ,♦ **N*»TAn 

MLft.HO-V Jf M^.r »G'Av- ft. *** • NUmOM 

Mfenort * CWst* 5Z 
IF IMSCVI - LVST*> 970*990*750 
9B0 1FATAJ. » | FAT Al. ♦. 1 
_ CALL rilPtl, IFLp* 

:'.A ’ GO* TO 99 $ U T * 7ff V 1 L ! MX ( B 1 » MP.OTl. LVST* 

99,0 NP V THX I t, 1 ■ 1 LV^TS • M 5 F V 1 ♦ 11 / 

' 1 F * ^EVTM X ( » 1 . L T * 01 ,„NC V T MX- C I ? ’ 

C 0 >»PUT£ BEAO iNOtCpS F 6 * TPAJMO 

XI I 1 n N Y rpG 
XI ?-» * Mr>tFi , | 

K| 31 « MCMAL 1 
K 1 ft I - MALflNl - I 
US > * Malgh 


■ V. 


M SI » MA| GUI <vi 
X C ■ .6 i ft MCML<*f ft 1 
M 71 • MCJ*L p l 




c 1 


K< n 1 • MCM1G1 • 1 

II H * M C »! T <? I 

XI 10) • MLpnt • I 

Hill * MLPpI 

KII23 * ( Mir,ni 

uni • n,l fin l 

inn • mc*fj» 'J; I 


vt'l;, 


(Y' 


-li.U'. 


■AV 


■' .•'••’I. • .■ 1 " : 


■\u. 


XI is’ ■ MC<WfJ'I i.'v '■ 5 : ' 

XI IA) ■ MCSP* # 2 • NFS ( | I • 1 

X C 1 7 ) - v* ? ■ * 

■ XI 1 B 1 « PC 5 m 1 ♦ > * mC« ( 7 ) « | 
XI 191 * mmcomi I i, "r’A * 




* v € € i ■> n^r r j 

X ( 23 IftMVtAB , 

X 129 J »HP f Art -1 
XI 25 »ftM« 1 ^iAn 
K 12 *>ftM*TAR -t 

X C 2 7 » ft'Mr'7 A* 
X123?PM#AB <*1 

XI 29 ) ftft«fT AH 



■ :i ■ r : 1 l' • 


!• J • orn-, • } 1 

' v \ y >< ;; ,• 

cr'i* * — 

•it .A, 



‘ i -i 


'I' , ) ' 

' ■■ ' '1 '( 

’ (Slv* * 


I \ *. f 1 

.'j. 'liA:’ !{.V« 


, h 


I k* . 








' ivrr- 


>i I' 


: ■ l. / Til 't 1 ' 

“Cn wrrf'i f'r-V 1 ...-' *.'• 

A 1 


?JS5» MApft 
Cl 55? I15|* 

'P,| 55;3 1 1*7# 


P r. 0 i- I I I * . 


ft. ^L-- ,-S^XAi ^.' .. .Av&ksinr. 


X I 3** »•*«♦*€ A«n *1 
X I 31 I ftML A MOT i 
Xl 3 ?»ft>i 5 UNP -| 

. - >y ». 1 1 wjmvwn . j.,-* 

■ r ■ ' 


ijiSMr : 




uull i 





M 


•••i 


•NPff 
• NC* 




v fl i / 1 ' ) a\ i j\ 

■' i. ' 


^ ■ .A':- 




.« , i '/,/■ 

"f/ * " «L'ftAftA‘ ft. a* •< l*>. *<>*»* 

! / - .1 ■» 


■ U\: . .u a*: . 




•'? ■WWW' 

• ** I M«M# 


~ » -•#^» '• • ^ ■ #<*■■ y «» ■»?<« . "■» ■« •»» •» ». ^ Jf*^**t* 

Table B-2 - Partial Lts*1j»g Suforoutine OAUXiV , 


OOIOC 

t If* 

CO 

CO IOC 

I2flp 

■ CO 

00 ion 

!?!* 

CO . 

CO 100 

I2?« 

CO 

Onj OP 

123* 

.. .. CO . ... 

oo I oo 

1 2<i* 

CD 

COICC 

l?S« 

CO 

C*> l 00 

I2A* 

CO 

***\ | 

127* 

. cn 

co iro 

12 B • 

CO 

Co 1 OM 

129* 

-o 

coinn 

1 3«« 

CO 

0** | np 

13 * « 

. CO . .. 

ecirr 

13?# 

CO 

OCIOP 

133# 

CO 

t r-n 

l 3«# 

CD 

coi no 

1 3S* 

CO 

roioo 

i 3 a# 

CD 

Co no 

I 3 t# 

cn 

CCIOC 

13B# 

CO 

CM J00 

i if* 

. . .CD 

:?no 

1 3n« 

CO 

C«|n** 

I3|# 

cn 

f cino 

ID* 

Cn 

COIOO 

I N 3# 

CD 

Como 

ms* 

CD 

ccsoo 

IRS# 

CD* ••• • .< 

oomi 

ms* 

SUBROu 

CMOS 

IP7* 

TMPLtC 

corns 

ISO# 

C.#;bF< 5|R ST, 

Coi OR 

1 S9# 

OOUfll.E 


_ „ ' «hcn' thc :. Oaux* entry isv^i 

A NO YPP IN ysTe A«^ 

POTENT T At * 'M*S S RATIOS* A#'; 
ON EVERY StlaS^flWCNT CAtt #£ 
_ POjNT oAVA IS CALLED, Tfl# ** 
that are Rf*n : Pieo are , 


j ory# ip» 

JfSp» #/ A^ r VMC^O|NTr»fS TON tNC . 
; *fefONS To V'np Ar-T PET UP* ■ ■• 
rti PHA.Sse CHANGE* -Yhe cntht • 
*atE4,E#AttQNS ON THE SPACECRAFT 


NpoF. CO ARE'/ - '• 

PERTURBING P'«cir E,5 OTHER TH*N CENTRAL fcOOY <S°BR BOOT 
IG<? (BURNS' {S'lpA* f RPPRN ) ' 

URP sunns TSAAWf LOPRR&.) ./, 

ATMOSPHERIC PtsT STANCE (SUB* DRAG) . ' •' 1M : 

POINT MASS EErrCT Of CENTRAL BOOT t COMPUTED IN OA|)XI 

n _ , . .£■ a * UAii.UAUA/. Cliflia C fc 


ROOT, THC SPACECRAFT POSITION IS ROTATED FROM M& AN»( 
AObE^T^THE^PNHAT IN f)AUX ARC THE VARIATIONS DUE TO THf 

central hoov. j? or The earth* ano the coao^h TtRM* 


• »**•***•«*•»**•■< 
unary* 1ST iYPP*VTAB,BTAe,ETAI»,i»TAB» 

ON 1 A*H • 0“Z > • INTEGER (I-N> ' 

SUBROUTINE OAUX 10/0 S//2 

♦BFV t CALPH6 *COAO?M,CCTUT #C6H« 

iltcDA02H( 2>»CGHR (121 
*». 2S> 

.PVHAT 


DOUBLE PRECISION *ASC *BFV ,C ALPHG »CDAO? 

OOUni.E PRECISION Cl»E »DAT 

DIMENSION flF* c 31 .CO AOjH ( 2> *C6MR 

0 1 HENS I ON UPT C *». 2S> 

DOUBLE PRECISION PAOH tPVMAT 

INTEGER PER3oO*RFLHRO . 

L»|HENSI0N NcENTP( 2)»NCS I 2I*NJ 

0 ,pE«?'tnni hi 

DlHfMSlON PvMAT I 3* 6) 

DOUBLE PRECISION SnANfiC* TB1.0CK .TOPERTiTOURN 
DPUHl E PRECISION TTOMt 

integer ST0C*0*VECJ ,TES %1 _ , 

DIMENSION TbLOCEI * 1 * tbperT i ?i»tbur 

f» * T P 0 T I *l,VECJ < 21 

n * ti.lth 1 3 ) . .. 

DOUBLE PRECISION RnPC. 3,31 
DOUBLE PRECISION TCCfF(3> 

OOiinLE PRECISION zi.1, 101 «TIhC2U> 

DIMENSION 1 1 IHC ? & I 

double precision TABqUTI ». 121 

DOUBLE PRECISION ROEC 1(3*31 . - 

DIMENSION TTOH ( 3* 3I,TT0HT ( 3* 

OOUBLE PRECISION COnST. CONFIX* SCRCOH* EBUF 
COMMON/ VENX/VOTAfl IS ,5*31 


f 2» »PA0" I 21 
»TORA« »TPOt ,TTOH 


TBURN I STYOPAS 


common/ /Const i2$o>» konst C 7S0>» 

c KONEIX ( SOM), length I 

e , KPOINT I 40). SCRCOH I. 

DOUBLE PRECISION THAjD 


CONFIX (ROOI 

100) » IPOINT I *0* 

200) » EBUF. ( 1000 



* 


0*?2S3 

? A 4 


Table B-2 - Partial Listing Subroutine DAUX (Continued) 



eeno 

235* 

1 

ceiio 

23A* ■- 

- 

eoi«o 

237* 

’ J 

eoiso 

?3n* 


5?|l| 

239* 

*■ 

roiB2 

2HO* 

1 

erstj 

?H| • 


E0IR1 

2H2* 


COIRS 

2-n* 


C?I ««6 

?1m* 

' ■ 

C91 17 

2*»e* 

1 

00ISO 

216* • 


Ccisi 

2 * 1 7 • 


CCI57 

ZHft* 


r?is3 

2 H ?0 

1 

f ?I5R 

25r * 


mss 

2S1 • 

’ 

CCI56 

25? • 


C? 157 

25 )• 

1 

C?I57 

25^* 


CO 1 65 

255* 


C? 167 

256* 


eci7i 

257* 

t 

roi7*» 

25»* 


C0I77 

2S9* 


C?7pn 

260* 


CC?13 

26 | • 

' 1 

cc?n«« 

26?* 

i 

CO 205 

263* 


CC?OA 

' 261* 

; . 09 

**" ?r>7 

265* 

( 

c??»r 

266* 

ro 

C?2I 1 

267* 

ro 

C0717 

26«* 


C?2|R 

269* 

1 

0*2 15 

27n* 


0C2 16 

27|* 


or? 1 7 

27?* 

i 

ro??e 

273* 

1 

CC??2 

77** 


C n ? 2 3 

275* 


' C**?7R 

27a* ~ 

> 1 

CC725 

277* 

i i 

?C? 26 

27b* 


CC277 

279* 

t * 

CC230 

?«"* 

► 

rc?3f 

2 A 1 a 

r ^ 

CC23? 

262* 


0C23S 

20 3* 


C«?36 

2‘Jb* 

- 

CC2R 1 

285* 

f* 

P32R2 

26 A* 


C02R5 

287* 


PC2R6 

2BB* 

' 

CC2R7 

289* 

. 1 

09250 

290* 

’ 

CC25I 

2 9 | * 


CC25? 

292* - 


ITTOHT Ml 

‘vecj i i i 


CONFIX ( 251) ) 

I * ‘VECJ Ml i LENGTH ( SO) ( 

I • ‘YES , . KONST ( 1 51 1 I 

C***EN«> STANOARD fOMH 0 N F 0 R SUBROUTINE D*U* 1 0/01/72 

DOUBLE PRECISION SUNaR? 

DIMENSION SljNAMr C 12. IS?) 

OOUBt E PRECISION Ma 5 maT< 3 . 3 >. ROMAT» 3 # 3 »; rOTmATO .31 
DOUBLE PRECISION Hill YPM). YPP (II 
double precision Trad 
double PHFCISIOn TOLD 

DOUBLE PRECISION UVKMAtIA.A) .TUL^VlMil 

DOUBLE precision «TAn<5 l M.VTAB(5,l)»BTA8M2»ll»€TA8<1,|l 
OOUBLE PRECISION CENVEh<3» 

DIMENSION TQAOCl) , 

0 IPEN 5 I 0 N BOI.pYO) 

DOIIBI c PRECISION EFMAt 1 3 • 3 1 f TEHPI3I 

INTEm OPT ( 2 I 


21>* 

? 9 <.a 


EQUIVALENCE ISRaNGE.TRi 

EFmaTC7»3j»EFMATI3»|I»EPMAT*3»2I /1*0#0«»/» 


/ 3 


18 ♦ 2 


J1 


I 


DATA EFMAT I | . 3) . 

| EFMAT ( 3 . 3 ) /| .0«/ 

Data opt/i?,i2R/ 

IFIN7TAB*EQ.0» GO Tq S56 
DO 555 J J« t . NVTAB 
DO S5S I 1 * | . 5 

. VOTAPI II »JJi"VTaB< I I . JJ ) 
e*£ CONTINUE 
554 CONTINUE 

NCOt. • < N Y« 3 1 
IALT ■ 1 
' I PH I NT « 0 
IOPWT ■ 0 
II « (KVEH.fl 
J I » 0 

‘ IF IFLDll I .p.STnCRD) . nE* *1 
KFUST ■ 3 • j i ♦ 2 ? 

II ■ IABSKKVEH.il*?) 

J) ■ 0 

IPAOM . FL|>I l I . ? .RPLWRD I 
IF (IPADH ,nF. a) J | ■ I 
I POTE ■ NFUST ♦' 3» Jl 

IPDtm * IPOTE ♦ 3 • VECjMl 

IRONY a IPOTM ♦ 3 • VECJ‘2l 

NFU5T • I Ron V ♦ 3 • (NbrSOL 

IEPSTR ■ N f 1 1 S T ♦ 3*(NSftLEP ♦ 

I E TOT a NSOlPP ♦ NCNSEP 

INARCH - IEPSTR ♦ 3»IET0T 

icc<*i 

ItRlTEM.Sl) IHAc;CN 
e? FORMAT I I OH |MAScN • 13) 

WRITEIA.S7) MHCON 

57 FORMAT ( 9H nmCON - 13) 

ARITCIA. 55 ) NY 
55 FORMAT! 6)1 >,Y ■ t3l 

XI - 6 • IkvFH . | ) ♦ ! 

LDPIGS ■ 3 • MaPTIa.KI) ♦•|APT|1»IJl ♦ t| ♦ I 
TOLD • TBLoCXM) ♦ 1000. DO 
TINJT • TRioCKM) 

CALL FLlP<NVTAB+3, IFLPI 


k nmRconI 

NCNOKPI . 


. RPITE(K0UT,fP7) M V T AB 1 I ,.|l . ! ■ | .51 f J»I.NvTAB | 

F *J . ,JLLy..l ; L. Vj* ?. ,* ?S* i n V r n< MM h<: • / f \ 5 c? 5 • I R I I 


S«02|4?0 


• NEW 

• NCR 

• HER 

• NCR 

• NCR 

• NCR 


• NEW 

• NCR 
•NCR 

• NCR 
•NCR 

• NCR 

• NCR 


C‘2-8 


Table B-2 - Partial Llstlno Subroutine DAUX (Continued) 


NOR* 

‘•O'** - 

Sift# 

** * » • 

«l ?• c 

in* c 

9M* C 

US* 

*»U# 

MT* 

«Ufl» 

«t|9« 

*i 2 t • 

H? 7 m 

N2*i» 

H?T 

H 7 t>» 

H?V 

M3"* 

**?!• 

*•33# 

*us# j 

15 a* 

h: 7« r 

4 »3s>* C 

•» 3v« C 

‘■■ir* r 

0!ARN OS t, c* 


HASmATiI tj» 
PVMATIf.JI ■ 
ISA * IS* ♦ 
CONTINUE 


• e.oo. 

PVmAT<|,J> 


rotate back to MEAN 


if t95o 

•Nf. til 


ROHATtltU) 


IP (NCENTR<kVEH| .NE, t I I 60 TO 275 
DO 273 IM.3 
773 TFMp (I) * TPOT ( f ) 

60 TO 2«« 

275 EFMaT < f .2 ) • EEmAT<?,|) 

EPMATI2.II m -ErMAT<2,j) 

CALL MMPY| 3 ,3,1 EFPAT,3 .TpOT»|.TEHP,| 
7«5 CALL MMPY C3.3, I , TT0M,3,TrPP, I ,TPOT. t I 
IP tuGPOT *rQ« ? »nR • I FLAG .LT* Ot 


Oft 286 1“ , 

CALL mhpY I 3,3, | ;pPMAT » 3 . VPPtNJCS’ , t . TEMP . | ) 

CALL MHPY(j J,| TT0H,3,TEMP,I .YPPINjCSI , I I 
286 NJCS • NJCS *3 
GO TO 290 

7«* CALL ROTAT C -OPT ( TROTpLI , TIME .TEMP , TPOT. <CTUT> 

IP I SOPOT .FQ. O .OR* I FLAG . LT • 0 I GO TO 2090 
DO 289 l*t,nGP0T 

CALL ROTATI-OPTMROTPL) .TlMf, YPP(NjCS). YPP t N JC S > . CEJUT I 

289 NJCS • NJCS *3 v 

.•9»* CONTINUE 
??9S CONTINUE 
79ft CONTINUE 

CALCULATE ACCELERATIONS DUE to LOW LEVEL THRUSTING 
ANr*/OR DRAG IF REQUESTED 

IPIKVEH.EQ.2t G 0 TO ?9| 

TNC. JCST TOR EQUALITY BETWEEN NON-INTEGERS may not be MEANINGFUL' 

IE« IVENT.EQ.3.AND.CftAD2M( 1) .EQ.Ot GO TO 3SQ 

Recompute attitude computations for 2 h or b»q mooes " 

0 N |. V IP TIME has changed, always recompute for #1 OR LVLh 
LOOK up VEHICLE mass ano CONFIGURATION. 

292 Call vmass<wtab> 

293 CONTINUE 

COMPUTE body ATTITUDE to ec| ip matrix required 




TO 290 


IFt lBATT.EQ.ft> GO Tp 291 
CALL BOOATT ( STAfl ,FTAB . Y . Vp I 

compute low level thrust acC* ip Thrust required 

I F « | VENT. EQ. D GO TO 291 
CALL SHVEnTIVTAb) 

ifinmcon.eq.ci go to 297 * 

IF < TDATT.EQ.ft. A hO.NmcON.GT.0) &0 TO 299 
CALL BODATTiPTAb.ETaB.Y. YP> 

CALL V il A S S < T A B ) 

IElNMC0N.GT.C-. AND. KVEH.EQ.lt CALL MASACClTPPl IMASCN) »MAS«AT> 


GO T ft ,->97 


IGINAL PAGE IS 
POOR QUALITY 


Table B-2 - Partial Listing Subroutine DAUX (Continued) 


C-SSP 

* Ac * 

C5SSS 

baa* 

CCSS7 

9A?* 

0CSS7 


ccsao 

*A9* 

C?SAO 

B7o* 

ccsao 

* 7 I • 

f -sap 

*7** 

t 6 i 

TO I ALNfl 

CfS6 l 

*7<|* 

f?SA3 

* 7** 

f?SA3 

*7e* 

ccsas 

97a* 

r?SA3 

H7ym 

2C r »f* 

*7P* 

rrsAs 

* 79* 

c?saa 

*8-. 

ECSAA 

*8 | • 

?“CSAA 

*87* 

CC5 A A 

**B1« 

fr^47 

*8** 

CC571 

*0S* 

r?S77 

* 8 * • 

r?S7? 

*87* 

“?S7? 

* 8 • • 

C?S 7 ? 

* B ♦> * 

CCS73 

••vr* 

f?57* 

*9 |* 

CCS7S 

*»*2* 

CPS7A 

•0* agno 

c?s?a 

* *s »* 

»*?S7A 

***• 

re* 0 *' 

HVc • 

CC*OI 

**A» 


*V7* 


**«• 

C? A** 

9*0*. 

C?AOA 


CO*OA 

SOI* 

fCAOA 

507* 

CC60A 

SOT* 

ccaio 

SC** 

CCA t 2 

SOS* 

??a i s 

SOA* 

C?A 1 * 

S07* 

CCAI7 

S**r* 

??41 7 

sc** 

£*?A| 7 

Sin* 

C?A17 

SI i* 

CC41 7 

SI7* 

CCA 1 7 

SM* 

CCA 1 7 

SI** 

C?A?0 

SIS* 

PCA73 

s t A* 

CCA7* 

S 1 7* 

“"A?* 

SIB* 

677 

SIP* 

CCA30 

S2D* 


GO TO 2*7 
7** WRITEU.290J 
2** FORmATISAH CORRECTION 
t * 

2*7 CONTINUE 


VENTS REQUESTED AND NO ATTITUDE MATRIX INPUT 


• NEW 

• NEW 

• Mf* 

• NEW 

• NE* 


CEuTRAl body ■ EARTHT drag requested? 

THr TEST pnp ECUALfTY RETWfEN nON-iNTEGERS HaT NOT BE 
2*1 If <C0AO2 H(kVEH» •EO, p.00 aORa NCENTR I K V£H I »NE* 31 
CALL DRAG C v» TP, SUNARY, JSZ ) 

SUM UP ALL PERTURBATIONS 

3S? CONTINUE 

KBUPN1 ■ KqiiRN ♦ I 
GO TO |3*0,3A?,37C> 




GO 


KBURni 


GET LOP 
3*0 


BURN INFORMATION 


If » IFLAG- fo*?! 60 TO 390 
CALL LOPrtRuJ YPPjLOPlSSn 
GO TO 3*0 

GET ICS BURN INFORMATION 
370 CONTINUE 

CALL IGSBRri YPPILOPIGSI ) 

39? CONTINUE 

THE test for foUALtTY BETwfEN NON-INTEGERS nay not be meaningful * 
IF IPAOMIKVFHI ,EQa 0.00 ■ AND. ( IPA0N*NS0LKP*NCNDKP) «EQ» 01 
• r. 0 TP 3*1 
CALL SOI «Ad I 
39| CONTINUE 

00 *?0 I ■ I , 3 

TCCfFCU • I-YI T »/TR3l*cMU*CGNRf ICEnTRI 
*?C YPP ( I I » l-Y( I >/TR3 )*CmU«CGMR( ICENTRI ♦ 

I AG < I I ♦ TBURNIII ♦ TRAC I I ) ♦ TLLTH<I> l< 


icEntr, kfust»nfust ,y,trao#Ypp I 


ITBPERTl I ) ♦ TPOTfll ♦ tor 


IEI lOFLAG.po. 1 1 GO TO *50 
GO TO *70 

ssr continue 

00 *5 1 I - 1,3 

CENVEHIII . TAB 0 uT(|,IO) - 


Till 
I l.OW 


• NCR 

• ••I 


INTO A UVW’FRAme rased UPdN L iu 5 L P 0 SITI 0 N JrT 

AS 


if att i turf more is si, tilth 
ROTATED into A UVW frame RASEf 
vehiclf amp vfmiclf velocity vector ano printed 
TLLUVJ. IF Moor IS not si, TLLUVW IS 5TAN0AR0 
liVw COORDINATES COp TLLTH 

*SI I E I MOE • N£ , * I rCNVEH(l) “ Till . . 

CALt UV«I C**UVCH . TP , UVWMAT I 

CALL MMPY(3,3,| f UVWMAT,A,TLLTH,l »TLLUV«,| | 

WRITE IFOUT,IO?*» 

fPO? roMHAT (73HrCKM60E PRINT FROM DAUX I 

write <KOUT,V*0J TBLoCkI I » .TPLOCRIRI »T0L0SV,ICALL»HDr,I8ATT # I vent 


rc**i 

C?AN* 

t s 7 


57,* 

S779 

f .?n 


99R P0RMATII3H tBLOckII, . p?&a|B,|3H TBLOCMW, . 02A.t»,|pH TOlOSV 
t026alfl/9H ICALL ■ 12, 7 h MDC • |7»9H IBATT * 12, ?M IVENT ■ 12 » 
ff R I T F I K n t J T .0*7) ( TPprJJJ » TBpTRT ( I) . TOR A C, I I ) . TBUR N I t ) « T L L T « < JJjL_ 


r 


ro 


cn 


PCMI 

C7691 

re**? 

C?*59 

rexs* 

c?*«s 

C?A*5 

CCA45 

ctaas 

cc*a* 

cr6&7 

r?*7o 

C*-A7» 

^47? 

re* 7? 

CC472 

r^c* 
rr7os 
ec7 |0 
.'’71 I 
??7 )* 
C**7|3 

rrzn 

fC7»7 

CT7JI 

r?722 

r?7>s 

PC 726 

••"777 

C?730 

P?7 3 1 . 

?C737 

Cp7 33 

PC*}* 

fp73S 
T?73A 
Cr.7*»(J 
7«»i 
CC7*»? 
rr7*i3 
cr75? 
PC7Af» 
P^ 7*n 
CC76I 
CC74I 
CC74I 
PC 74 1 
C?7*l 
C-741 
PC 7 4 | 
C?7*7 
PP773 
CC779 
ClOnS 
r:cn3 
C 1 003 


C10C3 


S? t • 
S 22 *. 
S 23 * 
S?M 
S?c* 
S?A* 
S?7« 

s;n* 

S?9* 
S3n« 
S3 )• 
557* 
S33* 
S3«« 
53*4 
S 34* 
S37« 
S38» 
S39« 
S*in« 
Si i « 
S“?» 
S«»9» 
5*1 * 1 * 

S*)t,* 

S**49 
5**7« 
5*t r* 
S*t*>« 

5 *» i • 
SS 7 * 
SS3. 
5S«« 

ss«;» 

SS A* 
SS74 
S.S4* 
s!s*>« 

54"* 
S 6 | • 
S47» 
$&?• 
5449 
5 6 

5449 
5479 
Sf 89 
549* 
S7i>* 
S7|. 
577* 
573* 
57*49 
5759 
574* 
577* 
S7fl* 


Table B - 2 - Partial Listing Subroutine DAUX (Continued) 


990 FORMAT! I 3H tBLOCRM, . 024*|B,t3H TbLOCK 1 8 * y„&2£* !SiA9 H . T ?fe°» * 
,026.18/9M f A ALL • I*,7 h HOC » I2*9H 1 0 ATT * 12, 9N IVCMT ■ |2.» 
J *RITE)K0UT,997) <TPOT< I ) »TBPERT) 1 ) ,TDRAG< 1 ) ,TBURNi I ) iTLI,™! | M 

997 V 0 RMAT ) I X4v|TPOT72X4HfePf RT2PX5HT 0 RAG*I XSMTBURN21X5MTI.LTH / 

‘•RITFlKOUrt !?*I ) * < TRAn J n tTCCFP< ! » «TLLUV«f( I J » |«! *3 » 

FORMATI IX.frKAO* ,2IX,»TCCFF* »22X»*TL<LUVW* / l3()2A*lfl» ) 


1*01 


e 

c 


CFMATI|, 2 » • -ErMATIl , 7 ) 

EFP AT 12,11 - •EFMAT t 2 » 1 ) 

CALL MMPT ( 3 , 3 , 3 ,EemaT, 3 .TT 0 MT, 3 ,RNP »3 
EFMAT»|, 2 » • -EfHAT( 1 , 7 | . 

EFMA T C 2 • t 3 . -<rrMAT( 2 ( n 


1*03 


»RlTE<K0UT,|9**3l (( RNP) ! .jl ,J-1.3) ♦J*l»3> 

FORMATIIX,*’ ftNP IN STANOARn FORM# / ) 3026 • I B I 


|«A 5 

999 


1 1 SI 


1 IF 0 


WPITFIKPUT, IP^S) ANCLE 
FORMAT (I X7HANGLC ■tl»*«|B/l 
FORMAT) 7 h yYPYpp/ ) 3030* I 8 > ) 

xjo » b ase ♦ tmlockhi»oav 

CALI .»ULCAt.(TfHF7»XjO) 

00 | 1 S| I * 1 t *> 

ITIPFU) » ▼ 1 ME 7 ) I • ♦ .COIDO 
SEC ■ TIME/)*) 

00 1 186 J«» .3 

3 . 403 *T«PFRT< J> 
3.4»>T*T8URN( J| 
3.a03*TCCFF( J) 

3. A03*TORAr, ( J ) 
3,403 *tLLTH( J) 
3.t03*TP0T( j) 

3 . a 03 »tpAD( j) 

4. f>| *yP<J> 

' 3.4D3*YPP)j) 


1191 


2C J,?> 

7) J,3> 

H J»**> 

Z ( J . 5 > 

71 J,4I 
21 J, 7) 

7 « J,»> 

Z C J ,9) 

71 J. 

CONTINUE 
7) I . I ) ■ T|MFZ< A) 

7)2. * ) • T 81 ncK C*0 /A. pi 

7)3*1) * Angle *57. ?9S779S| 3?823*f9r>2 - 

*R ME) KOUT . I 191) C I T I ME ) J ) * .)» I * S » »5£C 
WRITE CKOUT, 1192) 2 

EOrHAt ( • C nonnnCHFCXMOOE FrOM OAyX Af BMj 




• ••I 


S° 


t * 


• nredDnocnnonnnn* */) 


1 192 FORMAT) • 
hG 2 *> * 1 8 / 


SEC * • *F4.3*9X» , STFP < Hr ) ■ 


*.S()X.I2» *)X»F6.3 
• »G2S* 18»9X**«MA <0E6» 


2 

Ip 


TflPERT 

TCCFF 


tilth 
i * TRAO 
.» TPP 


) ) 
(Ew/H 8»»?1 
Cpp/HR**2) 
CER/HR.** I 
(fp/HR** 2) 


• »,3G?7.18/* 


* »,3G?7.»8/* 


f • 3 G 27 • i 8 / * 
* * 3G27 • 1 ft/ * 
» * 3G27 • I 3/ 


TBURN 

TORAG 

TPOT 

VP 


1ER/HR**2) 

1 ER/HR»* 2 ) 

<ER/HR«*2) 

IF.R/HR) 


♦ ,3627*18/ 
• «3G?7* »*/ 
» ,3627* )»/ 

• #3G27*I8/ 


JQ -0 

<=R' 

pa 


Rr 1 TE f ROUT » 1 1 9 4 j M rBEC f * L t » M«» » 3 ) » t“ 1 • 3 » 

FORMAT < *0 RbECI matrix in STANDARD f 0RM*//|RX,3G30*18) I 
ARlTElKOUT .1 193) l T ApoUT ) T-* 1 0 > _ • TA80yTI!»lll * 1***J ? 




119*) 

(193 FORMAT) IX,* SUN* ,22 X»*M00N» 


/ 1 2024* I 8 ) ) 


S7 9* 

r « «* • 














, 




I 

N 

• 

at 






01003 

S7 V * 

S8f»R 

C 


CIC03 

C 


0 1 cos 

58 | • 

c 


0 I 003 

587* 

c 

AOO 

Cl Of 3 

593* 

c 

fart 

01003 

58q* 

~ c 


01003 

585* 

c 

SET 

C 1C?3 

5° A* 

c 


- | r rn 
CIC?5 

i«: 


«7? 

C1CC7 

5*9* 



ricio 

CIOI t 

S9-»* 
59 | * 



?|CI« 

59?# 


4«n 

riou 

59 34 



ClC! 7 

S11* 



?1 070 

S9S* 



Cl C73 

59S* 


65- 

C 1 C 75 

697* 


44? 

01175 

S^'A* 

c 


?I0?7 

S»9* 



Cl 03 1 

40"* 

■ " ■ 1 


?i 037 

4**1 • 



C l 0 3 3 

607* 


. 

CI03R 

403* 



C103S 

604* 



? 1 0 3 4 

40*3# 



*10*1 

60** 



C 1 ?«*•• 

ACT* 


441 

P | ?«|R 

ACM* 

c 

C | 0M7 

*0> ACNncTIC* 

THF 

r ion7 

407* 


47? 

CICSl 

4 In* 


nos? 

411* 



C 1 OS 2 

617* 

c 


'“I'iSH 

4)3* 



C 1 0S«» 

41 «i* 

c 


rtoss 

4 1 S* 



ric',5 

414* 

c 


Cl OSS 

4(7* 

c 


?1CA? 

41 A* 


70" 

C 1 06 2 

419. 


75? 

61043 

6?n* 


C 1 04 ** 

42 1 • 


755 

ClOAS 

47?* 



P1047 

4 ? 3* 



C1C77 

4 ? 4 * 


74? 

C» 100 

47s* 



01 101 

424* 


??? 

Cl 102 

427* 



» ? •• 


• •»« M • ••<••»■•« IM» * ■ • # •••»» 


1 


Table B-2 - Partial Listing Subroutine DAUX (Continued) 


ANO THE C0A02H term 
SET CENTRAL »00V TErM 
II • IARSIkvEH-i » 


v ^ W >'VUI » fcfrn 

II • IA«tS<KvEH-| ) • IB 

IF IFLDU I ,2 # ST0CW0» .FR. fl .OR* IPLAG • 
TRS»TR3*TR2 

A ■ 3.00 • C M U * CGMr ( | CEnTR I / TR5 

A #l i. All i m I « 


LT« 0 ) 60 TO 7 S» 


- C M U 

do 4?0 1*1,3 
B II >*VC T > *A 
na t ftuTE 8 s i 


,V.t»PVHAT» 
IlCENTR) / 


IF t 
CALL 
CALL 
CALL 
CALL 
CALL 


1 FL * 
mmpt 
mtrn 

MM** Y 
MMPY 


*NF. 3) 
TEST 
> CO TO 


TR 3 


CO TO 
rOR 
470 


ire 


.f<l« Q » GO TO 47? „ . 

< 3t3,3,EFMAT.3»MA5HAT»3»R0TMAT,3 

< 3,3,EFMAT#3.MA$MAT,3 1 

< 3,3,3 ,ROTM#T ,3 • MASMAT ,3»ROMAT,3 
{ 3,3,3,TTOM,3»ROMAT,3,ROTHAT,3 ) 

. , i . - TT 0 mT» 3»ROMAT*3 ) 


J2 term 


MA$mAT,3 } 

,3 • M ASM AT', 3* ROM AT, 3 

{ l,3,3»TTOM,3,ROMAT,3,ROTHAT,3 ) 

MMPT ( 3,3,3,R0TMAT,3,TT0MT,3,ROMAT»3 I 
00 441 t • 1,3 
00 441 J * |,3 

PVMATI I ,J) * PVMATII.vl ♦ ROMATCI.JI 

TEST - ftC 


00 461 J * |,3 

PVHATC I ,J) * PVMATII.wl ♦ ROMATCI.JI 

test for cdad?n term 

TEST FOR foUALiTY peT*ePN NON-INTEGERS mat not 0C MeAMlNCfUL* 
IF CCDA02 MI fVFMI »Eq. n.OC> SO TO 750 

1 | * lKVC»-;-| » * |8 ♦ 2 * 

**• (FLD< 1 1 ,2,SToCf.'OI .EO* n * Tft 7Stt 


bc> so 
0) 


to 750 

60 TO T50 


IF 

1 fl as • (Flag 

00 700 1-1,3 


YPPI2I ♦ n • TnRAS ( 1 1 / C0A02M(KVEH) 
CONTINUE 

CALL HATPLY(PVMAT,Y|i»| # VPMI ,NCOL,VPr<RM 
CONTINUE 

SO TO 740 

»Y<!» ,YPC !» ,YPP( 1 I ,!•» »NT) 
♦ 1 
I 


IF< lOFLAG.Nr. I ) 
WRi TE f KOUT ,9*91 
{PRINT • IPmINT 
l O p R T ■ I OPR T ♦ 
RETURN 
ENO 


••■I 


rN° OF CONPfi *T|an: 

07*X ‘SyMrOlIC .. 

oaUx code relocatable 


5 D, AGNOSTICS. 


02 SEP 7» 
02 SEP 7« 


*>*|09I | I 
04109111 


n 02R3S406 I ** 

I 021/3112 « 

0 023 73*172 |1 


9 FOR,. MASACC.MASACC 

ii. .*'* 1 rn£rR*fJ V rvFr ! 1 I 


F V Fi 75 A »<F«Fr« I rVFI f l 73119101 1 


01 DEr 7 • 




149 


- •* • 


Table B-3 - Partial Listing Subroutine ftASACC 


ceioo 

Sn# 

CO 

COIOO 

si# - 

•“ CO — 

reieo 

52* 

CD 

CJIOO 

SJ» 

CO 

colon 

5« • 

CD 

COIOO 

SS* 

CO - 

CCIOO 

S4* 

CD 

cc ion 

57* 

CD 

or i on 

5b* 

CO 

coioo 

59* 

CD ~ 

Cc 100 

*p* 

CO 

C?| no 
CC) CO 

* 1 • 
*>• 

CD 

co 

Como 

*3* 

co 

mno 

6<i. 

CD 

CC | nn 

*S* 

CO 

Cl I 00 

*>*>• 

CO 

moo 

h !• 

CO ■■■■ 

CCIOO 

*B« 

CO 

moo 

*?• 

CO 

ecico 

In* 

CD 

Co I no 

7 | • 

CO 

C"t nn 

!?• 

CO 

cei no 

73* 

CO 

j r*o 

7S« 

co 

ceioo 

7S» 

CO 

CO 1 00 

7 *• 

CD 

CCIOO 

77* 

CO 

CCIOO 

7s« 

CO 

PC ! 00 

79* 

CO — 

CCIOO 

BC« 

CO 

C** 00 

«)• 

CO 

CCIOO 

8?* 

CD* • 

CC10I 

83* 


mp3 

w»i* 


CCI03 

as* . 

C***B 

ecios 

8a* 

C? 1 05 

B 7 • 


rcio* 

RB* 


mo? 

89« 


m ic 

9"« 


com 

7!# 


CCII7 

7?» 


m i 2 

73* 


m i ? 

7 M • 


cc ! I 3 

9e» ‘ 


ecu S 

9f,m 


ecus 

9 7 • 


m i* 

9 •• 


COI 17 

9<J« 

- — -- — - - 

coi ?o 

100* 


CCI2I 

!«!• 


CC 1 72 

107* 


CO I 23 

1 C 3 • 


ret 2 *» 

109* 


m 2s 

IPO* 


CC) 2S 

10A» 


eei2s 

107* 



STOCWD. SOLVE/CONSfO CODE WORD ToR STATC/ORAG 
KVFH - CURRENT VEHICLE 


OUTPUT 

calling sequence 

part - starting location roR vector or mascon partI^ls 
Masqat- mascon matrix 
common 

T POT , ACCELERATIONS oUe to GRAVITATIONAL POTE w# 

GmmCON- total mass of mascons in lunak masses 

CARO 

NONE 

print 

none 

tape 

none 

REMARKS ANf) RESTRICTIONS 

none 

SUBROUTINES REQUIRED 

VMAG2. MATOpS, FLtP 

METHOD 

MASCONc; ARE TREATED AS EITHER POINT MASSES OR *S ‘WASHER 
~ SHAPED* BODIES. IN EITHER CASE THE MASCONS ARr ASSUMED To 

oe rigidly embedded in the moon 

SUBROUTINE MCINtT(M A SCON,MASCW,*LPHA.GAMMA*VTABl * 

COMMON/T I HRK / H r ReT,THrRFT»TTRuI 

■GIN STANDARD C0«M0N FOR SUBROUTINE MASACC lO''98/72 

DOUBLE PRECISION Br y * CMU ,CRCf» *G«LUnT.OMMCON»SCRAT .TPOT 

INTFBER STtJC'VO 

DIMENSION BpY ( 3 > « C R C B ( U> t SCRAT I 200UTP0T f 31 

DOUBLE PREr r S J On tr^jd 

COMMON / TRAJEC / TRAJb c 20?)* ITRA*> < SO) 

DOUBLE PREClSIOw CONST, CONFIX, SCRCOM. EBUF 
COMMON/ /CONST ( 2 SO I t KONST <2S«3>» CONFIX fSOOl 

i£»jv, I POINT ( *01 

*70) . EBUF IM>00» 


*8 


c » KONFIX ( 5 OM, I ENGTH. * 1 CC* 1 * {POINT ( *01 

C» K P 0 f r { &0>* SCRCOM 1 C70». EBUF HO(jO» 

DOUBLE PRECISION T Tf ST , T I * T 2 , S , T S * * 

DOUBLE PRECISION CNt 

DOUBLE PRECISION TTEhp< l>,VT4B(*»t) * RBEC I < 9) * W£ I GHT * f&LOCK ( | ) 

DOUBLE PRECISION Hf.RFT,THEWET,TTRUl .ttrue 

EQUI VALENCeiRBFCI < I ) .Tra.JDC ISS) I ' 

EQUIVALENCE («Flr,HT,TIIA.|0( IS7I » 

EQUIVALENCE ITDI OCK ( f ) iCONFIX ( I70> ) 

EQUI valfncf ( NMCcON * I.FNbTH ( 60 ) ) 

EQUj VALENCE(NMCsOL*LENr,TH|79M • 

EQUIVALfNCFiMVT*B,K0NSTC227n 
EQUIVALENCE IREy I t » , TRAJO I l» J 

» » (CMij , CONST ( IS) J 

I » ICRcB II) * CONST I 28) ) 


• NEW 

• NEW 

• NEW 

• NEW 

• NEW 

• NEW 

• NEW 

• NEW 

• NEW 

• NEW 


Cel 2S 


(GMLUNT 

IrMuffVJ 


CONFIX 
rnwr t X 


I 2S2) > 
C 2S3) ) 




• ii , • O '/ ,•»' • /'• 


Table B -3 - Partial Llstlnb Subroutine MASACC (Continued) 


• :;vV v,;, 


3St*S 
231*5 
WIH 
3*125 
*3125 
m 25 
*3125 
33125 
3*125 
331 >5 
3 312'* 
3312? 
32122 
30127 
'M77 
*2127 

**| 7 7 

331 1* 
*3151 

32132 
“2 I 37 
32117 
3313.' 

32133 
32133 
32133 
*2133 
2*133 
*3133 
23133 
231 13 

* 2 1 1 3 
*2133 
2213** 
2*135 
*3135 
2* t *0 
*2P«7 
*21 •*? 
*2{5«* 

* 2 1 *• 5 

22M* 

2*1SC 

32153 

2*15* 

*2157 

3*140 

* 2 I 4 7 
*2143 
231*5 
2*1*5 
**!*(• 
3*1*7 
*2 I ft' 
2 * t 7 1 
2? 1 75 

;*?n? 

53203 


!?«• 

I ? * • 
I7t» 
174* 
125* 

1 3 ~ • 
131* 
I 3 7 • 
131* 
1 3«t • 
135* 
134* 
137* 

1 2 r • 
139* 
I *'2* 

I 5 | • 
157* 
HI# 

1 S«4* 
15«* 
I «A- 
1 '.7* 
l«P* 
|59* 
15** 
15 t • 
15 7* 
ISm 
15** 
1 5 1 « 
154* 
1*7* 
15e* 
I5«* 
i*r* 
!*•• 
1*7* 
1 

1*H* 

1*5* 


JCHLUNt ■ , ■ . 
(GHmCON 
iTFLAfl 
:<VFH 

(NMCOOE 
(NMf.ON 
(SCRAT ( 1 > 
(Stncwo 
(TPoT f|) 


CONFIX 1 *5*| > 
CONPtX 1 *531 1 


fTRAJ 1 
ICftNFTX I 
length 1 
LENGTH < 
SCPCOH 1 
K0NF1X | 
TRAJO <' 


191 / 
* 5 ) ) 
7*1 1 

n > 

1 1 i 

si » 


' / '■. - — 


e«*»FND STANDARD fOMHDN FOR SUBROUTINE MASACC ld/0V?< 

DOUBLE PRECISION HAsmaT(3,31 . . • . 

pOUBtE PRECISION NARCoN, PART* NASS. PHI • L»»J* To* •HR 

, SPM I , C PH | » BL*Ht CLAM» PMCN. TR* TR2» . R3 

*»ARlo»flR*0lu^CT,POLN(51 1 » POl.NP ( « M [sPHE • C0?11ARi){i5 1 1 ^ J?u ^ofu? J .41 

* * . AT SUM » StlitWH » IrS A V 1 ft TS AV I I ARS* V2 , ATS A V2 *t>ACC ( 3 j .TERNS *®GH( 3 # J I 
..prill , «MTX<3), ALFmA. 6AMNA. sL» CU» CL* SB 

DlPENsVoN MARCONI * 1 »H«SCW« 1 1 ♦AL^HAt 1 1 »GAHPA;< 1 ) 

D 1 Hf US | ON PART ( 3,1) , P H C N I 3 .tuPC t 41 

• .SPAR ( 31 , R , V ( 3) »PMAT ( 3»3» »SVE5 1 5* 


• .SPAR ( 3) • R * V 

• , 1 OCR » T 11001 

EOOtVALENCf . (SrRA 

• , I STR AT l 7 1 * P H I > , < S C R A 

• ,( 5 CPAT 1 * 1 , 5 »*H t ) ! (SCRA 


• ,tbCPAT * .CLAM \ ,|5CRaT T»PfjCN 

• ».< SPRAT * ll».TP 7 1 , ( SCRAT < I 5 J»TK 3 
*,( 5 f«AT I I 6 I,TFMF ) .(SCRAT 


0 (SCRAT ilSll.lSCRAT) ,{SCRAT J l»#MSS 

1 , (SCRAT ( 3 ) *L AH ) *<SCRAT 1 5(,0 

) I (SCRAT ( *)»CPHt 1 » (5CRA,T 1 

\ , (SCRaT ( 9 1 .PMCN j * (SCRAT J 1*>#TR • 

1 ! i SCRAT I 151 » TK3 ) .(SCRAT I t5» »?COR 


Sl.lMC 
Hi ,HCST 


1 » (SCRAT ( |2» »TR 

.(SCRAT 1 51.GC0N 

> .(SCRAT i ZC'.RIY 

) , HSCRATl 3> »*00N0 

) . ( fSCRAT 1 
I . 1 1 SCRAT V 5 7 » I 


• .(SCRAT 1 21».PpAT 1 , (SCPAT ( 321.SVEC J . ,, * »0N A I 

"'*!» ISCRATl Il.tvPAR 1 , (fSCPATI 2J.1CJ7 1 »{!|CRATJ ?J»5°° N ° { 

• , ( | Sf R A T 1 4 >,rol>CS 1 ,( t SCRAT ( «* J » ! *i I 2 S 2 M 1 o! * l I 

*,( 1 sr«A 7 < 7 >,|i»AR I , t TSCR AT ( RtiMCST 1 .(ISCRAT* * J * I \ 

*. (ISCRATt |A),P| HTEs) / 

00UPLS phcmsion tihcoa 

EO'M VALENCF « COfiFI X C |70l .TTmEODI 

OATA PAXTR.i /72/ 

0A1 A wraRS / 9 / 

DATA HAS5AI. / .FALSE, f 
|N|TtAl.1 7 •- SURPOUTINE 

:cc«i . ;■ .- A ,v ; a ■ 

IF ( NMCON. GT. 2. AnO»NVT AO. C0«0» OqO*SORT ( TV y 1 
. 05 J»1 .NMCON , 

DC » I • I » f J V TAB 

. *>( j«t i *t ♦ a 4 

Jr « vtIm 1 1 *1 !gE,HASC0N( jT.YL> • BHD •V.TAO *;'#!* I > |PVU>|0 TO t 

GO TO S 

1 0 <: •*. t • I J-ll.SM 

lo :o.?o i j*i i *3** - ••••'/ •- 

IDt’M J* ( J-l | *9^3 

VTABI 3, I l«VTAB( t. 1 »*MAseOH( tOUMlT 

VTAB (5,1 )• VTAB(5» 1 1 + H ASCON UOlfH*! • - ” ~ ' 

VT AR 1 S » I ) • VT*8 ( 5 . 1 * *HASCON ( IPUR3) 

5 SST^SOl, (HASC»rCfAl.tA«*l.*IHCO0El 

** FORHATl VH HASCj. * 0121 * - '• 

NH*9*NHC0N 


2*705 

•~7.t 


RRfTFU .50*1 ( RASCONi t?T) • !**•> »MHl 

foR** AT < »5 h mA5£aN * 9,t>»l«.A.) 






» • #» .•.*♦•• •*•-■■ •» 


..»•* * »n«w * • ** 


Table B*3 - Partial listing Subroutine MASACC (Continued) 


»*»• 

UT» 

I AP» 
169 * 

! 7r>* 
t 7 t • 
!’?• 

I M* 

| 7 * 4 * 
lie* 
IT#.* 
177 * 
| 7 «» 

( *•• 

; • *• 

If I • 

I 07 * 
M3* 
l«*i* 
IRt* 

i «*• 

I ? 7 » 

J «f»* 

r<* v • 

i»«* 

i’i* 

197* 

I'M* 

■ I v«r 
! vn* 

IV/,* 

1 9 » • 
|9n« 
i9o* 
IS*)* 
231* 
237* 
203* 
2 r,t ** 
*?<;• 
23 '.* 

20 /» 
20 «* 
20v* 
210 * 


•mreiAfSo^) «HjseoHi tin » i«r*t »nh> 
ro*HATUOH mA3C0 n ■ 90t3»A> 

ENTH* MASACC^^T*"* 5 " 4 !,* _* AkaETIIKN 

ir<N*CCf'N»FO»e«ANO*NilCSOU*fO* 0 »HETUII ,l < . 

| If f l T ■ 0 

tn'2*HHcsot.*n M <cnN 
00 7^33 11*1 • 3 
00 2003 

WANT U 1 *jj>*e.i?oe 

T?£si»t)A 8 S< t«tOcK l*J»l 
DO ? r 0J JZr»l iNhCOOE 

K0DP0»M*SCw< -I?7 I ■ ■ 

rOONn-FLOl?* IP*K00€S) 

»r to*??* j*7 4 * * 2 *itooesi *ee.o> so T0 2000 

I XNT • 1 <Mt ♦ | 

I nur*»» <KoriMo*< * • *♦* 

ifjha 5 c«n( *6T,t2) 55 jjj 1222 

|M^ASC0N<!n'JM5j.tT,T I «0 TO 20C0 
IF IMASCOMC t DUN* I * 6 T*T 1 1 SO TO IT*" 

T?*r 1 

00 TO 1 99 | 

* TS»HASCOHUOUHhI 

1 ,!!;H?ssv!WS!i*riT>??W" t0 ”’ 1 

iso TO J 9 T 3 

» S.|T 7 - 7 SI/tTESt . 

3 CONTINUE 

CMT*.SS39S93TO-2 
9 |E/r-IKMT 


!0UmA*3* 1 JK7** > ♦! _ 

rA«!‘i»t§ZF>* 5 !‘ 52 fSi 


2 ? 0 i 


3300 


^>A«TI 2 ,|FZF» ,, S* tgocrt 

f‘A u 7(3»fE2C)*S** R8Ect< 

CONTINUE 

CO»*T I NuE 

retiikn 

UCTUHN 

ErlO 


O»;MA>f*Gf0HTl*CNT 

duha*i >/ wf 

0UH&*2> /wCiSMTMCNT 


•nee 

•NEW 

• NEW 
•«►%! . 

• NEW 

• NEW 

• NE* 

• NEW 

• NC!» 

• NEW 

• NtW 

• NEW 

• NEW 

•NEW 

• Hi* 

• NtW 

• new 

• NEW 

• NEW 

• NEW 

• NEW 

• NEW 

• MEW 

• NEW 

• NEW 

• NEW 

• NEW 

• NEW 

• NEW 

• NEW 

• NE" 

• NEW 

• NEW 

• Nf * 

• NEW 

• NEW 

• NtW 

• NEW 

• NEW 

• NEW 

• NEW 

• NEW 

• NEW 

-2SS 


SS 



f NO OF (OHWIlATlnNt 
iTiAEC C STHWOLfC 


27 0tA6N05T|CS* 


mU; £c code wccoCatasue 


» HAT TO 
28 HAT 7* 


win?; 


0 *W *«|76 

02011*70 

02011 * 3 * 


ii ,5 I JSR 

IN 13* 


01 OtC *W 


c ro*»* hcnpwc jHCwwwf ^ ^ t t , r „„, «ci iuvPf* * *vfi ri. 7fii00lfl» > 


/ 




si TOVd 



e*'o* 

0*403 

"cion 

ccio* 

*£907 

rr*ip 


c?“An 

rr^is 

f *HH 
rr4A5 
c«9 aa 

r”«*-y 

r.***»o 

?3*77 

** 47 ? 
Zt 4 * 7* 
CC472 
CC474 

*ns7(; 

r?v76 

rc '477 

ns**a 

ct.soi 

T*RD2 

f, r't.ry$ 

rcs-iM 
r?*07 
*: r b » i 
C?5l2 
** 51 3 
6^*s i -• 
;jsi5 


CPSI7 

C?52l 


??*?? 


3459 
34*9 
34 79 
3**9 
3499 
3**9 
?<»• 
3«>?9 

34 J* 
35-99 
35$* 
5s *9 

357* 
3Sn* 
3 S 9 • 

H"* 
34 | 9 
3* ?*■ 
343* 
36«« 
34s* 

3649 

3 4,7* 
34 a* 
3*?» 
37*9, 
37 | 9 
3719 


Table B-4 - Partial Listing Subroutine MCNPRC 


00117 

3239 

fin 3 

3219 


3759 

c*«U 

37A9 

ro*»»4 

- 3779 

TOUT 

37*9 

fo**?o 

3299 

ro 4 *?! 

SSa# 


31 | 9 

r*?i?3 

3 3?* 

CO^il 

337* 

CC*«33 

3399 

fOi3G . 

3Sr,9 

r^no 

33A# 

Z' A*»| 

3379 


3199 

con** 

339 4 


2 A»>9 

7 

3H|9 

CC AS2 

3 4 J 9 - 

o;iS3 

34 >9 

f?SM 

3499 


CALL T*»*lT <PT*I 

NUMSN m NUK*M ♦ I 

BRAD • "ASCONIS, JOASRl I 

Arab • hascqni 6,jwasri 

TH* • HASCoNf7^JWASR3 t *70t 

I PR » Q i 

CALL ALRHAI tf».THK*ARA0»*AA) 

CALL ALPHA 1 ! tPR.tWK*ARAD,»A*> 

I PR • 1 

CALL AL***a< tfP»THK» A»7AB»AARA J 

CALL ALPHA I | PR ,T HR » BRAD, BAP A I 

ipr • o „ ; ' ' . ■ 


CALL OAHHA(|PR,tMK, ARAD,A6A| 
CALL 6AMHA<tPR,THK,RPAP t B6A» 


CALL GAH«A< jPR *-fHK »BPAO ,BGA 
JPK > I ' ' ... 


CALL GAM«A< iPR *THK * ARAB, AGPA I 
CALI. GANNA ( f PR .tHK ,BrAD,B6PA> 


WALL li ( j?-r • Tn«, ,n*AU • nur«w 
*RITP(RnRUM|) ('AtPH(L) ,L“l ,3*J> 

CONT« HUE „ ; 

If INHCIO ,r<9* e> ®0 TO 340 

DO 350 1*1 ,nMC \ - — - 

17 • ROOM > 

IF: I 60 TO 340 

IF. 1 I? • <9F * I .AND. 1 2 »LC * NHC > 60 TO 310 
IfRR ■ TERR ♦ I 
CALL FLIP! I flFLPl 


CALL FLIP’ I # IFLP.I 

FCAnATlRO^rRRCR ILlE®AL VAL^E ON HOP <ARO VAtOE«»,It0) 

♦ BP Ml +'rl * t- #9 . 1 1 ■ ■ V- ’ - .. 


G6 TO 3S9 

cc> n t i our 

00 3*10 K* 1 .9 , ■ . 

IF IMH5C0N? 19RK *\ 1*1 .f«* P> «0 TO 3*0 


I PACK. - 2 
GO TO 330 


KASrOHlK <1 ?» • .mASNOmIr , 
k«con • Max (NMCOM, i 7) 

TRACK' * I , v •: V- * ;. •: 


NRScdnI I 4*k i-\ T>) 

CONTINUE :•, 


** I R ACK 


CONTINUE 


CHECK |F 5/C fiRCATjr* TmAN NKCON 


NhCdOF 
NMCPLM 
NMC50L 
IIHCCUN 
IN IT • 
LOC * 

1 • e 


CALL DECODE ( 1 • . m » TP A» I , 

IF <1PA*> V70,R?C3375 
-IF INI I ) •Gr .MNCnN? tJp TO ‘>05 •• •'.•••• •••• 

IbTt^p l R0& ( t N I T f - 1 » ♦ B » 32 > * t AO'S t *• f ***( I PART *21 

K • 1(1 I ) 

!PARI>! A8S( iPAR) ' • • " ■• ■■■■ - - 

IF‘NMSCONl|if9|PAR«*K>#EO«C> .60 TO 3?S 




IF I I N IT • E 0 • Mill! 60 TO 37S 
I N ! T • N< 11 

\°S ; to' * t 

».:;..rnnr . * i*mfdoE +1 


' 4 


}•» r t/4 ¥ "l» t. +* 

''t*. I,.. <. r A f .\ i. ijL ..A-:.:’. 


■ -}V ' ■ • . A , .. 


■ ;i > 

f i : h i * ■ 


• <• '•i'W* *' v V V ' k V ‘ V»« V‘ # a ' • i*i' i^" 


eesi* 

CCS21 

* 0^23 

CCS2V 

e?S 2 S 

r*5?A 

«i»r.?7 

'“r'SJO 

r??» 3 .v 

">eVl 

cecis 

CCS 34 
?•' *> ? 7 

rcsvo 



eesv i 

- • 3 A A • 


rftSR? 

3»7* 


rn^t, q/ 

3vr* 


f H 

3V»* 


“eSH-i 

372* 



3 V j • 


T’ISPS 

3»«i* 


e?5V4 

3Vc« 



3 V a* 


"•kjc, -» 

3 V 7* 


?'-5S2 

3Vn« 



379# 



• Op* 


z r s r , h 

VOt * 

OB 

r o 5 *» >» 

vo?* 

• 




ressv 

1099 

ro 

Ce655 

VCf.* 


C0557 

VOA« 


r?s4i 

**C 7* 


ee?47 

in«« 


A«l 

V' 1 ?* 


r«s A *i 

11 "# 


eei.' , .5 

111 * 


*• p 5 7 r» 

H | 7* 


0?57| 

VI 3« 


res?? 

111 * 


coaoo 

Vis* 


C04HO 

V I 4* 


eeAoo 

VI T« 


CCAOP 

Hta* 


ce/.o? 

V I ?• 


fe4cr» 

V ?«• 


•‘ew i v 

V?t* 


??t 22 

V 2 7 • 


e-A3i 

V23« 


C«4VV 

17** 


?04VS 

V7$* 


rc4s a 

V 2 A • 


er aso 

V 7 7 • 


C045I 

V 7 A * 


C3oA2 

V 79# 


CS463 

• 3n* 


-y«-v J . ,V ; ;: :-•■ ;; ■>.;,; 

:/ V^;,, .;■■■ 1 Y': : ;.0' ■ //'.' 'htity. ;V- r V 


■;?/'' r-r^trrr 


Table EM - Partial Listing Subroutine MCNPRC (Continued) ,■>" 


.'j ; 

,,V:' 


IF I IN I T • EO*N ( | J > 60 TO 375 
I N I T « N<n 
LOC ■ LOC ♦ J 
1 " 0 

npcode • nhc^oe ♦ i 

CALL PAC*<* t t'»,M5CO0(NMCOOE)*NM»J 
CALL PACA<79»7 iHSCOD(NMCO0EI • LOC I 
375 I * 1*1 

IF IIPAPI 3flV.Vn0.3fl5 
JR? J * 7 

nmCccn « nmccon ♦ I 

60 TO 3 VO 

3«5 J * I - 

NMCSOL ■ NMfSOL ♦ I 
3V* IBIT • 10 ♦ 2 • JABSUPARI - 2 

CALL PACK I IPfT.»»M5COOrNMCoOE> »J> 

IF ttP«P«NE»l> NMCPLM • NMcPLH ♦ 1 ... 

60 TO 345 


•k- : y\&- ■; 

>: • . \\ ■ ■ , « ' . 


ERROR 


3V5 IENP ■ IERr ♦ I 
CALL TLIPI f ,IFLd 
WRITE TxOlfTlaCilJVl N<il 


30 1C FORMAT I 73m fPRlIR ' - CANNOT SOLVE OR CONSIDER THE PARAMETER SPECIFl 
• EO FOR MASCfiN NO* * TV* t2H - HOT lNPUT> 


CALL DELCO'o 
GO TO 345 



r>fV OF HASCOHS 

TEST F ff « CwECKMODE PRINT 


VO? IF (MMCON.f/>,C> 60 TO VQ?> 

if ;,i list it y;r.si go to • see Sh 

calj fl»pir*?flp> ■ -V:., ■ 1 1 ; ■ " 

WRITE IK0UT t 4 »9'*?> . ^ i ! 

«" 5 ?.ro«K*T <75 h CfCMdpE PRINT FROM, HCNPRC^ , VX.I0HHAS3 RATIO* *8X 
, nHj.ATIT^^FtlVx*’HL0Nr,lTuOE*19A»6HALTlTUOE I, 1 i , 
VJ- 00 V2C J»*. ,*|PCON , : ’.:r. ; | >'/••’ <-rV’ 

CALI, F|. 104 | , IFLP) ji ! -A- 

ttRifrncouT,.rte> (mascqni i ,j» * *7> 

FORHATI7027, !0> - ;| , 1 1 ‘ ti-n;-, , ■' ; 

12? CONT t NUC . ; ■. 




WR JT€ MASCONS ON DRUM 


50** CALL RaNLasIKORiiMI tNoRoMl^I > > " v 

white ( kormm t i ( i.mascon it »!»,t ♦npaRSI , j«i »nmconj 

APITF<KO«Um MMSCOOIl J t | | , NHCODE » 

WR|TF(XO«Unt « ( (NMSCOnI {•J)« f a t5.2t) »J*I *NMCONI 
WRITE! KOUT ,vOl I l f ImASCONC l » J > , I " t . 7 > • J« 1 ,NmCONI 
4««« | | FORMAT ! 7f> t « . I ~ > 

call flipimhcou/7*h, jflp) 

WRITF IKOUT.FOe?) 

K »lflMCOM- » 1 / V ♦ l * 

WRITE OCOUT.Mei^I MhMSCONI ).t I , J-16.2H .1-1 .K) 

. . CALI. FLIPInmCOOe /7 *■ 2.IFI.PI 

WRITE (K0UT,H0V*)1 NMCOlJE » < HSCOO ( U » I * > • NHCOdE » 


.. •NEW 

\ VNEW 
\ •NEW 

.. \**o 


FORMAT(10H i|M*?CoN » /( 7(2* t 0l2)l> 


I * (M/MI - 

V - ■ • 


a 

w 


eetoo 

®A* 

e?ioo *" 

S 7 * 

re ioo 

SR* 

r"ioo 

59 * 

moo 

An* 

moo 

«l* 

moo 

A 9 * 

mno 

. AY* 

como 

As* 

rr ion 

AS* 

crloo 

A A* 

ec ion 

A?* 

re i 

A •• 

**—«/“. e 

A*f* 

men 

*•*« 

r**ieo 

7 a* 

Colon 

7 ?* 

r ? 1 no 

7 Y* 

moo 

7 s* 

moo 

7 s* 

moo 

7 A* 

?"HT 

77 * 

rp t no 

73 * 

Ef i r « 

77 . 

**C 1 00 

R n * 

re ieo 

S|* 

r oino 

«?* 

C* |oo 

BY* 

meo 

8 s* 

—M |f,n 

« 1 * 

r« i no 

« %• 

men 

6 7 * 

men 

RB • 

moo 

« 9 * 

moo 

BO* 

moo 

9 f* 

moo 

97 * 

r - 1 nn 

9 Y* 

*• I r»o 

9 s* 

**C ! no 

9 S • 

e- roe 

96 * 

ceioo 

97 # 

mm 

90 • 

eel oo 

9 R * 

so ion . 

ICO* 

«*- t on 

10 1 • 

CMcn 

107 * 

mt »0 

1 0 Y* 

C**inp 

10 s* 

cctei 

ies* 

mc 3 

n** 

mo* 

107 * 

cjios 

i?«* 

erics 

ICY* 

m*A 

I to* 

mn? 

III* 

m to 

117 * 

SCI 1 1 

111 * 


•■'ll? 


1 I • • 


Table_B-5 - Partial Listing Subroutine TRAJ 


co 

cn 

co 

co 

co 

er> 

co 

CO 

Cn 

co 

co 

co 

CO 

CO 

co 

CO 

CO 

cn 

co 

co 

cn 

co 

co 

cn 

c» 

co 

cn 

co 

co 

ro 

co 

CO 

1*0 

CO 

co 

CO 

co 

co 

cn 

co 

cn 

co 

cn 

co 

co 

co 

co 

co 

cn. 


CARO 


NTT*S - NUMBER or TIME TRIGGERS 

NYT«G . NUMBER or OEPEnOENT VARIABLE TRIGGERS 

TblOCK- time block EOR TRAJ 


TAPE 


none 


none 


coerrlcigNT 


OUTPUT 

c * LL NoSe sE * ueNcC 
common 
CoaP 
i ns vf . 
i*sv£ - 
in ag - 

jpnT - 

kfi_ag - . . ... 

TBlOCK- time BLOCK roR TRAJ 

CARO 

_ none 
print 
none 

tape 

none 


• OR AG . 
STAB POINTER 
WTab POtNTER 
CONTROL elag 
etab pointer 
control flag 


...ITTaLIZEO HERE 

INITIALIZED mere 
for INTEGRATOR 
INITIALIZED here 

FOR integrator 


FOR 

FOR 


BOOAYT 

vmass 


tor bodatt 


remarks 


ANO RESTRICTIONS 
NonE 


subroutines E 9 i i I R E 0 

TRJOUT (TRJ-.VRT) , 0 AttX , ENOSTP • TRlGER 


— J 1 ETM 00 


Traj is a 
solution 

EQUATIONS 
ING AN El 
TWICE* A 
obtain in 
select jon 
halving a 

ERROR. C 

OF INOfPE 
ultimatfl 
of differ 

ISaAC-NEW 
AT SPECIE 


double precision numerical integrator for The 

F A SYSTEM OF SECOND-ORDER ORDINARY 0 1 FFERLNT I *l 
• IT IS BASED ON THE COWeLL-ADAMS METHOD, EMPLOY- 

gmth-ordfr backwards difference formula sunMto 

B I I tl P a M I i V • a • _ A ■ I n * LJ — A n - n (? V . n TPB ■ «• 1 A * — - * 


OF 


non 


KNDtN 1 VANIABLt.5 ARE n|T» T K A — 

Y PRODUCES A TRaJFCTORY TAPE CONTAINING A MESH 
ENCE LINES WHICH MAY LATER BE PROCESSED BY IH£ 
Yon cOMPINATION roR I NTrRPOLAT I ON/ i ntegrat ion 
IC VALUES OF THE INDEPENDENT VARIABLE* 


c*« 


SUBROUTINE TRAJiTTRf,, VtRG, Y,YP,YPPtOlF»SrVENT,VTAB) 

nouoiE precision vtabis.u 

COMMON/ VENX/VOTaMS, 50) 

IMPLICIT DOUBLE PRECISION I A-H ,0-Z » » I NTEGER <I*N ) 
•BEGIN STANDARD COmMON TOR SUBROUTINE TRAJ 
OOUnLE PRECISION CpLOcK » TR bOCK 
p REC t S 1 Om COao^MlFI .COAR 
PRECISION CONST, CONFIX, SCRCOM, eruf 


to'os/fr 


• NCR 
•NCR 

• NCR 
••-I 


double 

OOURL E 


w u v i nr. r t c I J w ri ». ur '3i| vu«r i k. r. rup 

DIMENSION CBLOCKI S>»NTT*<; < S?f,TBLOCK( 


1) 


COMMON/ 


/const <2$e), konst confik hooi 

■ * * ■ .. .. . ... *.#.•« % ■ l • A* • 




_I»ble_B-5 - Partial Listing Subroutine TRAJ (Continued) 


oe*t* 
CO! I* 
0?*l* 

cons 

771 19 
771 is 
*?!|* 
CO! 14 
0 0 » I * 
-£l 14 
7-1 f6 
00! |6 
C01I& 
CO! I* 
-0116 
-Oil* 
00116 
0-11* 



7-1 1 A 

; 33* 


-7U7 

139* 


7-120 

135* 


-- 1 20 

1 36* 


C?1 71 

1 37* 


--!?? 

1 ?9* 


--I 72 

ii** 


C-S 7? 

1 8o« 


cri 23 

19, p 


--•23 

1 92* 


— » 73 

Ml* 

€0 

77 : 25 

1 1«» 

1 

771 ?s 

1 9^* 

ui 

- f I 2 e - 

Mfc* 

ro 

on 27 

197* 


r-i 77 

199# 


CO 1 27 

199* 


-71 31 

1 Sn* 

< 

-7131 

Mi* 


**7:3i 

157* 


*7133 

1 5 1* 


rci 33 

I5«* 


-- ! S3 

1 55* 


7-1 33 

15a* 


771 39 

15?* 

. 

77139 

ISA* 


--t 39 

» 59* 


77 1 39 

1 6-* 


771 >5 

1 6 1 • 


C 7 « 37 

162* 

V 

771 90 

16 3* 


77191 

169* 


--1 99 

1 65* 


-'*196 

166* 


£7150 

167* 


77151 

16a* 


C715? 

1 49* 


77153 

17?* 


57159 

171* 


.cohmon/ /Const < 2 so> 

C* KONptX « 

C » KPOlNT I 

DOUBLE PNeciSfOt, TRajO 

Double precision tilth cai 

COHHO* / TRA JEC / TPAJO < 
EQUIVALENCE tCBiOCKtl) • 


/const <*so>. konst <2so>* confix (9oo, 

XONptx ( 50l». LENGTH < 1001 » I POINT 1 *0» 
*S2 rN I * . SCRCO* * *00 > » EBUr (1000) 


UFtAG 

< CO A R 

< 1 vfnt 

1 CD A0 2H ( | » 
I JP‘lT 
t 1BSVE 
t ! *sve 
CKFl AG 
«K VFH 

tNTfNG (|> 

< NT * * 

< M V T R 6 

< TSt.OCK t I ! 


zrj'ji, i 

CONST 
T TRAJ 
TRAJO 
KONST 
CONST 
jtraj 
ITRaj 

ITRAJ 

ttraj 
KONF f X 
LENGTH 
length 
length 


ITRAsI 


1 ■ .. *T8t.0CKtl> • CONFIX { 170! ! 



^T f 1 > *VPPt| ».Ofr* I3*| » .TTRC.I2* t ! , TRr,T«2! i 

c iiS*^oAr?{2ie5^' ,0 ' -ss «& a«< v ». s . cs .*». v > 

DATA «SS!Gf 1 J *1»1 •»> / 9751299,00* 10-81956,00. I 0957!»29 . 00 1 

- I J08?00'J-.Do« 11391670*07* l»T7Si)R0.00» 1 2690^20.00* 

2 1 33 »5 677 • 0" , I 3 3 -5 4 "" • 0-/ * •' 

DATA <GC,*M( I > ♦!;*«•?> / 2782753. DO* 2190f»39.00* *207VNC.0O» 

1 2399070*0"* 2530080.00* 2721*00*00* 

2_ 3??R?CC.Dnt 362B0i7.no/ 

Data <STGS( ( > . f" I .?> / -337157.00. •37A068.00* -932976.00. 

! -5«163S.OO, -58.399C.DP* •665280.00* -665280.00* 

2 7»Do, | 33056"C.D0/ 

OATA -IGAMSi | } ,!•« ,9>/ -57281 ,00* -67906.00* -0750". 00* 

i '1 -z££' , * 0 C* -136-10.00* -191520.00* •3o2‘*Od.OO » 

-3628800.00/ ' 


data cssigm > , i«i »v> / 9751299,09* ic 

1 icsooo-i-.oo* 1 1 39 1 6?0*00* it 

2 1 33»5677«0". 133-56"-. 0-/ 

DATA (GGAni i ) , I 4 1 ,V> / 2-H2753.D0* 7 

f 2291899,0*1, 2399070,0"* 2 

2_ • . . 3-29577. on, 3628097*00/ 

Data Istgsc t > , fat ,9> / -337157.09. • 

I -5-16-9.00, -583990.00* • 

..... 7 o»t».u,h S( , '”sniT:ISi. 

1 - 1 7356", oo» - 1 3671C • 00 » - 

? - 67982 -. 00 , -3628800.00/ 

LOGICAL #FLAG ,FFLfcG 1 

? W I T I AL T 2AT t ON PROCEDURE 

CALL TRJOUTfOfF) 

INITIALIZE pRAG, VEN t, AnO ATTITUDE HOOELS 

«F< IVENT.EQ.il GO TO N"0 
. CALL SHVNTl * 

9"7 CONTINUE 

DO 925 1*1,3 
925 JLLThi 1 )*0*o" 

I F 1 Co A D?H ( | ) .Lf» 0> GO TO 603 .. 
C0AR«CDA02filKVEM» 

6"1 CONTINUE , 

J P N T * - 
. 13SVE-1 

1 KSVE* | "■ 


10-81956,00* 
I 19750*0.00* 

2190039,00* 

2530080.00* 

•376068,00* 

-6652«0.00» 

-67906.00* 

-191520.00* 


10957529.00* 
12690-AZ0.OQ* 
^ * 

220779", 00 > 
272(600.00* 

• 9 * 

-932976.00. 

-665280.00* 

•02500.00* 

•302900* oo»; 


77IS5 17?* 



ST3LK9 • TniOCK(9) 
tf r r c ▼ m . 
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CD 

t 

in 


u 


run 

s?m 

0*156 
0*157 
0*1 At 
C*H> 
r ?l A** 


??ii 5 


P*|7fi 
rri?i 
**i ?? 

?* 1 7* 

eo 1 74 


**• I 77 


"?? 0 I 

?"?r? 

r*“?o j 

7 n* 

-•' 7 ns 

cr?o 7 


17 ?* 

17S* 
|7«J* 
175* 
17 a* 
177* 
1,7 n* 

1 7 ?* 

I?** 

bk- 


i*5» 

1 Pa* 

t 67 * . 
I ft n * 

1 *•>♦ 

1 V«* 
I’t* 
IV?* 
I 9 t* 
IV** 

1 »a* 


*02*5 

*“'' 7/7 


2 3*. 
7 V j , 


t#T?. ?-5 * Partial listing Subroutine TRAJ (Continued) 


? 

c 


i;iir : i BL0tic ‘^ 


ENO or INITIALIZATION FO* ORAG.vENT. ANQ a TT1T(?0E MODELS 


IF (KFlA6«FQ< 


ir iT-VS^TT^St^' •r R * NTTf? 6 (KvEH».L€*|l 60 to 1000 
lZ-ril r * C ll * 7 1 *t T * OoOOl 60 TO 1000 
* ^T-R^C?i2i * f 0i flCX 111 
II 5. kjtxg / Tmlock! i )> 


■>» , 


r*r?i«5 

1 V'A* 

rr ?» ? 

197* 

rc?ts 

i 9r,* 

*7?lA 

199* 

tT7 | 7 

?0n# 

c???n 

• 70 1* 7 

Pn?a| 

2C** 


201* 

r*'7? , i 

20** 


7 r 'G* 

rr7?*» 

72/* # 

*r77A 

207* 

* , r?ti 

?C»* 

r*'?.72 

209* 

# ’*‘?33 

21** ' 

C?73G' 

711* 

*?737 

717* 


21 T* 

P**? M | 

2 1 H • 


715* 

CC?*3 

214* 


21,7* , " 

TT7*S 

21«* 

P*'?*A 

219* 

r??*? 

270* 

r?7S0 

22 | * 

C?7G3 

227* 

c??ss 

223* 

cr7s* 

22*f* 

r*“r57 

22 V* ~ 


22a* 

CT7AI 

727* 

r?7A7 

7? A* ■ • 

C07AH 

229* 


• ^(#61 *#* 7 ' n . ■ ’ • MIvva \ I | j ■ // /• 

ir;, S7tu v : Y ^p^dl T ?i!c*i ? f* ?? ^ Tjs ^ T ^^ veNTJ 

'i ■ ; = 

70? CON V I Mi/E : ‘ ■ '■■ ■•: ■>■'/}:.■■ '■ 

KFLAG-3 < 

■ ■•; 1 flag * ♦ j 

I^LOCKfVl-T^LOCKfVl/D.ftOO r, ‘< 

f Jr? rTWSisft I ■ tTs 2 (S?TI 7 

i CONT I NUf ' ’ 

XFLA<j"0 ■ ' '• 

? 5 "„* 


■o 


.. 4 U‘- 


fNDSTF 


KC « I 

r,r u » • 




«ll> 


TRGTf 
• NgAR * 3 

It ‘I'YitlO) .NFftA - NV 


, I VI » I J ' 

Do LI TO 1*1 .NT, 

01 Fi!,t) • yfFIrJ 


OIF <2# I >*VI I ) 

>lt « F « * * * - 


111? DIF » 3,1 9«V p If ? ■ 

IF ICBLOCKi?) .»fe 
.... «*« * NENR ** 

K 5 • 3 
I ! Ft AC ■ * | 

5AVH « THLOcF f A) 

ASSIGN 1 13* TO (fOK 
I »?* OP A B 5AVM/?#90«t 
T6L0CK C * ) *$A VH * 

KCNT * I 
60 TO 7C0p 
MS? 00 I I AO I*j ,NY 

oirJ 5 fn , S , ;fJ,«° HW 

GO TO i | ft?* 

| » «»*• C I F | 5 a 1 ) ■ v I | ) 

US* *«!’ " 0*Fl?.I) 

T P I | ) m 0 I F 13 » ] i 
1 165 TPPC | ) * 0 1 F ( 1 1 1 ) 
TBLOCA(I) ■ TOLp 


.!■ i. 


• fBLoCM f l > 60 TO 1230 


■Vl., 7 ,' 


/m •* S - 






J^^CNT.EQ.?) 60 To 1170 


"V'-.ll'l * * 4V. . V. 


- 1 .6.4 1 . **>7^*r **•*>*•<• <- - - »• — ■•. 






; 7 ' ■. .: 

AaL'i i .rj L» > ill’;. Af. -' 


.♦HEW 


•••I 


Table B-5 - Partial Listing Subroutine Continued) 


OB 

I 

VI 


cotas 

cctsa 

C»7S* 

f"7sn 

cotsi 

C0762 

€0753 

55755 

cars* 

ta**l 

r? 7*5 
0*749 

C"7 VS 
CVAA 
fi*77i) 
C~77 I 
C*7 73 
**7774 

r iooi 

n-oj 

*• i - *« *; 

~1“14 
r »oo7 
C 1010 

cm » 
0101 2 
rmi 
? i a 1 4 

ri o» 7 

oio?n 

C 1 071 
«* 1 077 
01073 

C 1 C?9 

r I not 

^ J w l ” 

*1374 
C 1374 
ClOJO 
>1050 
Cl 130 
? 1031 
01037 
CI03S 
01017 
5100 
0|O*iO 
0|0**7 
0 1 C *• 3 
, ClOO 
' 010*14 

0»0*«7 
CIOSO 
C10S0 
OIOS* 


IS** 
040* ‘ 
*14 I * 
96?* 
943* 
944* 
965« 
04 44 
*147 • 
9AV* 
069* 
*»7o» 
*57 t « 

9* j* 

V7 3# 
•57*** * 
97^* 

*5 74* 
17?* 

9 7 a * 
979* 
9fl«*a 
94 I • 

*»e ?• 
**- 1* 
9# 4* 
985* 
984* 
*I37» 

. *5 6 r. * 

* Hi’* 
9V«a 
Hi* 
M** ?• 
S»l* 
9 9 M* 


a ** 38 TBtOCKtU • T BAR 
# ’ NOS • NY . 

ASS ION 6559 JO K0§_ 

.ASS l GN 4000 TO KEXT 

60 TO B0C9 „ \ 
in** KFLA5 • KSavE ... 

TRBT'MI " TfttOC«<ii> 

IF l*Fi. AG*LT*2'1 60 TO 6050 
vcrcp ■ 8 




•v~: 


,:V ; ' 


if nr fi...-,. 

KST£P 8 

TByOCKMl 


4****7 


w ■ t IM_<- - ■ 
tblocxT*»» 

NYSAV2«NY 
M Y » fl Y S A V 1 

ASSIGN 4 095 TO k wT 

60 TO 9000 

4 **I 1 S IF !KrtA6*E0. 31 60 Tr) 4fl*57 
60 TO 6096 _ . « 

1 F 1 mVTa«*F«J.C 1 60 To An6 1 
00 6598 Jj^ifNVTAB 
DO A? *18 U*li5 
VTA«I I!,JJ)-V0 TaB11I # JJ> 

4-«8 C0N3!Nijf 

A**A 1 CCr-TlMOF 
RETURN 

4**94 Tfl'.0r.K<9)«m 

»jv«ny-sav7 
GO TO |**W0 

***«;♦» ir ( kflag»eo.**> «0 To 4«s7 
* " 1 F 1 1 • Y T R 6 •LF- 5160 To AQS 6 
I F TR » I 

ASSIGN 4055 TO PRESET. 

60 TO 6500 
4*55 1PTR - 2 • 

FFL.ir«o 

A 5 M 6 M 40S7 To K r ESET 

G n re, A soc 


• NEW 
•NEN 

• NCR 

• NEW 

• MEN 

• NCR 

• NER 

• Nf R 

• NER 

• HER 

••-2 


H ? C # F 9 J -• ■ -Y >1 '• ' - » • W r%- — * 

mi* Gn TC A60C 

seT;«».g 0 ?; , S5fj«» "* r " 0T e£ " t *"»"‘ r ‘’ L ‘ 

99 n# 5 4*57 ir 5 TptOC< ( | > •E9»TR6T<2> 1 6 ° TO ®01S 

499* 60 TO 6005 

50*\« C FNn TIME ONLY LOGIC 

soj# C begin scan r 

S0?« . A*AO SUGAR « “ 1 • f)0 

IF « YTR 6 < ill » Ilt*3»C0*J1 60 TO 60*5 

DI*?SmTI£. H.Y HOT Ot MI»N|N,rUt. 

•mar T,C IF IO.EO .0.MO51 60 TO 60*5 

S07i R « YTR615, I> -.YTRG<2*!» 

IF 1 R*D » 6 T • 1 60 TO 65*5 

I Ft ■’ 


S07» 

S09* 

SO*** 

S10» 
SI 1* 
S I ?• 
S13« 

S 1 9 • 


« ■ IT K'l » 3 ( I ' . 

IF 1R*D*6T»4»?D*i! 

YTR6<I*I> « !*0j 

s ■ n / <R-n> 

SU'MR • OMAXt (S,SUBaR) , 

A**AS CONTINUE 

c CNn |^ C (i^JJAR*LT•®•bP" , 60 T ° 60 IR 


10S9 

1 5 ?.S 


SI 5* 

51 A* 


ASSIGN aoa a to kkIoe 
GO TO AACO 


- n R * I « • •„ 


<•*>- • ' •* •• 


■ *«kW'W "• • " 


??J ?2 

fCJ*S 

0 ?i»* 

r?s?7 

CC 3 VJ 

0 C 231 


r*? 33 ? 

^0133 


o\ 


co j: 

— 335 
r?i j 7 
CO 3*10 
CO'>*M 
003*17 

— 343 

— ■45 

— i ** a 

003S7 

003*i3 

?C>55 

— 354 

*r 357 

— 3 
00341 

— 34 * 
“5 >44 
'■037ft 

— 37t 
00 '77 
'•0371 
r « '.70 

— 3 7*; 
03377 

— <*r, t 
r * •• r» | 
00 **o 5 
**0‘*00 
COM 7 

— 4|2 

C— I 3 
00** 1 0 
CO * I o t 
e— 1 4 

— a « 7 
C — 75 
0*»**7 r » 
r—?A 
00*27 
CO** 30 
C0**3 1 
—**32 
fci**3-i 
CO* 4 3 3 
CO 4 * 3 3 
P?*<35 
00*04 
CO** 37 


09 **«f» 


SC — I 


J»5* 

374* 

377* 

37»* 

379* 

3»«* 

341* 

3 n ? • 

3') 3* 
3««» 
3f*S* 
3BA* 
3»7» 
3 AS* 
;« .» o • 

34 A* 
39 | « 
39?» 
39 3* 
39*4* 

39 ** 

39 A* 
347* 
394 * 
399* 
A — * 
4C I* 

*tf*?* 

*403* 
4Q*I* 
4?5* 
4*JA* 
A 07 * 
AOA* 
AC*>* 
MO* 
91 I* 
*•17* 
413* 
41 *** 
4|$* 
<t|i* 
417* 

4 J A* 
4 | 4* 
A7' 1 * 
47 t* 
3?7* 
47}* 
42 a* 
47S* 


Table B -6 - Partial Listing. Subroutine TRIQER 


I1JG0 TO 3|0 


3*« 


3»*» 

375 


ispect • t - m 

NCOLO • NCFNTR(ifVEH) 

IF I NCFNTRIKVEH) ♦NC. 

NCCNTRIKVEhJ ■ 3 

YTRKl|»3» • 4.0*» - 

VTHf,(2,7 > * CSPmINI 3) 

TTKGI?,3> • CSPH I N I I I ) 

6fl TO 350 

IF (NCENTR(tfVFH) .Nf. 3)60 TO 3?0 

I r < M ,F0. *»G0 TO 325 

N C E N T R ( K V E H I • |» 

YTRr,« 1,3* * -I • rtO 
VTR6I ?,?l • CSPH I N ( | 1 ) 

GO TO 350 

IF* MCEHfR(rVFH) «Eo, I0>60 TO *40 
MC ENT R I F. VEH 1 ■ 10 

00 330 T*2,|0 - • - • - 

J* I- 1 
YTRkI I , T * * * 5*00 
33* Y7H*i I 2 * I I ■ CS^HlNIj) 

GO TO 35? 

34** J ■ 4-1 

NCEnTR (FV fhI ■ j 
YTKr,I7,7> » CSPHlNtj) 

00 345 1 ■ 3 , 1 0 

YTRG* 1 , ! > ■ -! «n? 

IF * J • ?E • 3* 60 TO 350 
VTRG* 1,3* ■ C*D~ 

TTRGI2.3* • CSPHtNIJl) . ... 

35*» COnTTNUF 

icfntr ■ ncfntrikvch* 

CCHRIIi* ■ G«L*Jt|T 

IF * ICFNTr .EO. It » (tO To 3010 - 

IF t ICENTr ,nE. J J GO TO 3090 
If ( K 1 6 • I. r • > f,0 TO 3090 

*F*D EARTH POTfNTJaL INTO VSTR 

CALL RaHSET~< KORUMl , N0RUHI7* > * 

RK*t> I KORUMl * < IV5T »(I* , I *. K 1 5_ • K ) 6 „ | iU „., u ..... 
CAUL GPP T a I IVsTRtHfWEJl ) , I VSTr < MC«E J I * 3 ) * I VSTR I Mc» e C 1 1 I 
• IVsTRifljEll • IVSTR<HCSCI> > 

GO TO 3?90 

3*»|*» if * K 1 8 • l T # < 1 7 * r,0 TO 3??0 
RE *0 MOOH POTENTIAL INTO VSTR 

CAUL RANSET < KORUMl , NCRUMIB) > 

RF AO ( rORUMl » < tVSTRII) , I * K 1 7 • Kt« » 

CALI GPCTa ( I VsTRlMc^MJ* > • I VSTr I mCRM J I * 3 ) 

I IVSTRIMJMI) , !VSTR«MC5H|> I 

3**7« CONTINUE 

JOflO CONTINUE - 

C6MRI I I » “GMtl'NT 
J«9« CONTINUE 


345 


• !V5TRtNc»»«CI » * 


• NEW 

lir 


• nE* 

• ME* 
••■a 


,0«44*^sTIC* 3 The Y TEST 7 foR ep'* ALrtr*PETSEEN C KN5iNTE6ERS MaT mot |C 

•177* IF ITMINT .FO, TBI.OCK < | * »ANO, I RCQ ( NC0I.0 1 , CO, 3 ‘.AND. ICEnTK 


••?»* 

47** 

430* 

431* 


|.fG, ITAMOO* 60 TO 375 
ITAROO • ICENTR 
1 RE6 * NCOLO * • 3 
TMINT • TBLOCKC | ) 


437* 

437* 


CALL rPHEHl 

CG H R * I I * * n»> UnT 




RAAiMtMHkRi 









k 


; / 


< 

Ot 


N 




C07bO 

fl!MI 

55797 

5?**9A 

55997 

55950 
C?«<4« 
??*• 45 
C5970 

?5«7| 

5*97 3 

5*977 

5*3 6 3I 

?c!3? 

55S0? 

55573 

5' , S?«* 

J757* 

rv*?r 

556 3>T 

rr^ji 

55133 

f_l37 

5'»135 

5511 * 

ecu/ 

"oS*» I 
55*97 
' , ?15t 
55551 

55*53 

5*559 

C*5e,* 

5*5,5? 

5*iai 

05*63 

e**3A*i 

??v.5 

C55A? 

55*71 

55*77 

5517s 

55S7A 

55*77 

5?' 52 
55*03 
55* |5 
C06 I A 
0*6 I 7 


5** ?5 


«?• 

733* 

735# 
1)U 
137* 
7 3n# 
737# 

Sin* 

•»•»*• 
77?* 
77 3* 
77«* 
Sic* 
79 A* 
77?# 

mi: 

os**# 

15?* 
•I5>* 
9S*i# 
75c# 
7**# 
7$?# 
7S9# 
7$## 
St«* 
7*1# 
76?* 
si>» 
769# 
*45# 
96 { • 
767# 

9«o* 

7*?# 
77«# 
7 7 1 • 
S7?» 
77?# 
••7m* 
775# 
A 7 A* 
77## 
979# 
*179# 
A Ilf*# 
*0 1 # 
7«?# 
*103# 
*•*<<• 
705# 
A«A» 
907* 
98 A# 
909# 


99"* 


Table B-6 - Partial Listing. Subroutine TRIGER (Continued) 


37 f 


370 

305 


CALL fPHCHf 
CGHR! I 1 1 ■ gHi UnT 

rVEHl ,E 

|f‘l.i>II'»HF*rESJ6» TO J78 
FORMAT! )AM*eCNTffAL BODY ■ I < > 
EPMjN a T BLOCK ( | > - TMtN<ICVrH> 

RR f TE ( 6 * I I wAMEfl I im.PomIm.TRI 
CONTINUE 
OO 38? t*»,3 
Till • YUl ♦ T*0OUT(I f MCOLO» 
VP!l> ■ V P < j ) ♦ TABOUTt !#3*NCOLD» 

«tot ■ msoi r p * nc«sep 

XMUl . c M U * cGmRMO 
ir I N T 0 T ,«r. 5) CALL 

“ “ - TO 


CALI fPHCKf 

CGMR(ll) ■ GMLUnT 

lFtMCENTM|«rvER) ,EQ« | ! 1 CGMRC I t 1 mGMLONT 
|F<LISTTT«MF»TES>CO TO 378 
FORMAT! » AH*» C FNT»?AL BODY ■ !<> 

— • 0 CKl| > - TMfN«KVrH) 

mAHE(L> *M,CPM!N.TBLOCK(Rt »!TUI,yPI r I ♦VPr<n l l«|,3> 


• NEW 
•••I 


305 


370 

37 1 

J»3 


{» 'VI 2 ' «mf* VI CALL FPHACClNTOT.fCENTR.NcOLO.t.rr.XHUM 
1 F(l!ST»T*!JE.VE 9»G0 TO 3B5 
CALL FLt p !?,tFLr» 

fpilfi i^V^LftCK! |f E - T ?M|NlKVEH) * 

CALL FLIP! |9* IFlP) 

JerLAR ‘ K ° UT » 1 * NAME <L» ,M»EPM!N.T0LOCR<7» ,(Vfl),rP!l) #YPP» 1 ) , 1*1 ,3) 

TF (TBLOCKI9) «GE» ?.f>(J) GO TO 370 
KFLAG ■ 7 
GO TO 371 

i i » nvfhsf!kn :h» 

! »LE. NFVTMXfrVEMl) GO 
r ' - nvehse»kveh> - 


IF (II 
NVEM5E! fVEH) 
Ge TO 3’S 


TO 

I 


373 


Oil I 'I 

Y * » IfIVEHSF IKVfH) #71 . | 
SCVf NT III) • M ♦ ISPECY 
IF* L ! 5 T I T .NT* YES ) GO 

r * ii ri 10 / • tn ni 


TO 379 


• • i- 1 j ■ • • • ■ • » 1 l i uu in 

CALI FLfPlj.lFLP) 

PR ITE<KOUT * f 7A7 ) J|, 5FVFNT I III* 

1 707 FORMAT! IM«J,73X»7IH*** SPECIAL EVENT 

* ft, *i | 8) 


377 


ISPECY, H » T0|.OCKI|» . 

NT / t X • 1 10*015*6*21 I0»5*» 


37S 

•too 


T0LOCK ( f ) 
YTRGl 1 ,18) 


GO 


-|,00 

TO 701 


701 


7*»7 
710 
7 1 5 


1 3 


026. 

ISPECY • 0 

SEVENTH!#!) ■ 

IF IM ,E0« |0) 

RETURN 

IF I LNOTR iKVCh) »F0. 0 ) 

NVFMSCIKVEh) * NVEHSElKVEH) - i 
GO TO 376 
ISPECY « 2 
JK • ? 

IF ITBLOCXfn) ,(_T • O.Or) JK • 1 

IFIRTVPF ,EQ. HRTYPE(2l »07. RTyPE *E0#HRTVpE(5»» GOTO 707 
IF IfriTFL .GT. I) GO TO 710 
OCLfMl jit,kvE«) » TBt.oCKlI) 

60 TO 715 

OCLTMI JK .KvrH) • TP| OC|C < l ) 

OTL IMI JK ,<VFH) ■ TBlOCkH) 

IITCNT • ITCNT ♦ 1 

if <1 istit .nf# rrs) go to 7I6 

CALI. FL I P ( 3 , 1 FLp ) . * 

WRITE !<0UT,!.T) OCL 1 m » dTlIM 

FORMAT IVH-nCLlM ■ * 70? 3 • ) 5/7H ©TLIM • ,7023#IS) 

EpM ! N • TBLOCK C 1 ) - TNINIKVEH) 

CALL FLIP! 10. IFtP) 





• R I Tf t * OUT , | | N A ME ! L 1 * M * EPM I N , TBLOCK I 7 ) , I Y ( I I • TP I I I ,Ypp f I > , f I , 3) 


. ./ . ... t ti. a* 
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Table B-7 - Partial Listing Subroutine TRAJRD 


GO 

I 


) 

) 

> 

) 

3 


05209 

c:»i 

c-?o? 

;??oj 

0520 * 

1 5 ? r . S 

r?70A 
C9710 
59711 
55?l I 
5*171 1 
P-71 I 
5? 7 | 7 

tV?n 

r -772 

C 5??7 
5??3P 

5 5?33 
-?733 
55733 
5-733 
f?7*ll 
C57<*3 
557**** 
--?«•* 
59?«R 
C-7R* 
C-757 
C5?S** 

m 753. 
-?7«iS 
C 5 ? *» r » 
557*3 
r-7*s 
557*4 
557** 
-5?4 6 
552*4 
C-77R 
5-77 A 
--777 
T5777 
55 77 7 
55777 
5-305 
553-S 
5— '*■>* 
C5»0S 
— 30* 
55 3 10 
583| I 


23** 
23 1 • 
237* 
233* 
23** 
23s* 
23** 
237* 
23 a* 
23?* 
2*9* 
2«»1* 
71J» 
2*7* 
2**m» 
2‘IS* 
?*U.* 
2**7* 
2*4* 
2R4* 
25n* 
25 l • 
257* 
253* 
25i* 
2SS* 
75*« 
257# 
25** 
257* 
745* 
2*1 • 
247» 
243* 
2*1* 
2*s* 
744* 
247* 
24 ?• 
24«« 
77-* 
27 | • 
277* 
?7y* 
27** 
27** 
27** 
277* 
27** 
277* 
2 An* 
2«f* 
2«7* 
2*3* 
2«f* 

ia«« 

2*4* 

287* 


C 

c 

c 


c 

c 

c 


c 

c 

c 


13*21 £ ._ 

" IR • *7 ' 

IS « 14 
14 - 23 

17 • 22 „ 

150 CONTINUE 

y TRIGGERS 

IF t K2 «Lr« El I Go TO 7C? , 

CALL RANSEt I K n RUHi , NOROmJIII 
REAO ( K0HUM7 > t VstR<!’ t I * El 

AEtGNHFNT CONSTRAINT TABLE 

»n- IE I ICS *LF. K3 » GO TO 300 

m i*sis u 5°niT*,?rv M l , sM 

ALIGNMENT ( *OLC 

jn- IF < K4 »LE. E5 I GO TO EM . . 

• CAM. RANSCt < E n RU«7 ♦ NORUm \ 1 3 1 » 
READ ( F0RUM2 I < VStRII* * I ■ E5 

eo* burn constraint table 

S»»— IF I E8 *Lr. K7 » GO To $09 

CALL RANSET t Kf)RUH2 » NORUMjfRl I 
RE*0 < K0RUM2 l°( VSTRIII » I - «7 


K2 I 


KR » 


K4 > 


K8 I 


C 

c 

c 


jGs BURN CONSTRAINT TABLE 

5«n IT C K 1 0 »LT * Er » GO TO 45* , 

57 ‘ CALL RAN5ET * EdRUM? , nORUMI 15)1 

REAO C KORijm? » 1 VstR(I* » 1 • Et » EIO I 

EO* burn table 

in« IF ( K 1 2 «LE • K « 1 I GO TO 70? 

* CALL RANSET I EDRU«2_, > N0RUHj IGl^j t 


C 

c 

c 


~ re aH V*6*M l u rvsTRH) 
tGS BURN table 

• M n ir ( k 1 1 .if i Et3 I GO TO 8n*7 
" CALL RANSET* E n RUN? i N0RUMJI7) J 
REAO I K0RUM2 I c VSTRM* • * * K I3 

potential 

„*n CONTINUE 

EARTH POTENTIAL 


E 1 2 I 


KIR I 


llu*. Ui s lV T 5 tilm s ? SU*??. 

READ I KORUN* 1 I VSTRlfl t I • 


) 

•CIS • E 1 4 ) 


55311 




• .ii anfruTt ii 


Table B-7 - Partial LI s tine. Subroutine TSAJRP (Continued) 


i 

* 

M 


£931 1 
(93»1 
9*31 1 
c?m 
99321 
9933* 
9-322 
99333 
99372 
9*319 

9*312 

99333 

*93"! 

-«l*.3 

9-351 

9-353 

99353 

*93*0 

£9350 

£-3*9 

fKjAn 

C*3i3 

993*3 

99333 
9-33.1 
99331 
99131 
9933 3 
99 **0 1 
99*03 

c-*o* 

99“53 
--M0 
9-« t A 
99*30 
9-*»7* 
99*33 
99*30 
99*21 
99*23 
99*2* 
99*2* 
99*2* 
C-H7* 
99*32 
C-*3* 
-9*15 
J*9**35 
99- .35 
99-35 
-•mnI 
99**5 
P-4*A 

99**5 


I«M 

2«*« 

!!?• 

iv,: 

»**• 

2*5* 

2 VA* 

in: 

299* 
39<** 
-39| • 
3??# 
):»» 
3?«* 
395* 
335* 
333# 

5 *** 

3*7* 

3 I -* 
311* 
31 ?• 
313* 
3 1 «• 
3t*» 
31 A* 
317* 
3 I **• 
3»«* 

32- * 
32 1 • 
32>* 
321* 
37*.* 
32c* 
32#.* 
337* 
32** 
329* 

33- * 
331* 
337* 
3 3>* 
31«* 
335* 
335* 
313* 
3 3a* 
3 

.1*9* 
3* i * 
3*7* 
3*3* 
3*** 
3*5* 


C 

C 

c 


*00* POTENTIAL 


•*« If I KID *LT» K»7 » GO TO I,??}. 
’ ' ' CALL RA**ET * Ko*W*l • NOWU^ljl 


C 

! 


READ ( KORUHl 
*A?S C0*5 


I I VSTRI!> t 1 ■ 


> 

*17 


» K|8 1 


tC 20 .L T • K»* > GO TO JlPO 
RAMSeV (KqRuMI «NnRl|M< 91 » 


••*9 If f 

* CALL 

RE AO 1 fORUMt » « V5TRII » » ! 
ir 1 K27 *LT • K7l » CO To tOSO 
READ ( KDR jm 1 ) < VSTRU* » ! * 

• •5- If U3 A*LT.k^5J |»0 TO MOO 

CALL RA w 5ET|30R»I*I »W0RUM(90I 1 
00 I960 I l-OS.*!*!*** 

1 1 I At “ W ♦ 6*2 

(-6* R£AO<KORUH|) (VSTRiDt ! * II 


K|7 i K29 
K2I ( K22 I 


I m- Continue 


S 


VTAP/BTAR/ETAR/RTAB tables 
vtao 


If <K7*«LT*E73) r,0 To 1 70$ 

CALL RANSE TCKORu«2»NoRIJMt7S» » 
REA0IK0RUM2) < VStR H I , t*K23 »K2* 1 

BTA* 

12*9 IE IIC26.tT.iC7S) fiO To I 10 ? 

1 CALL RAMSCT(*0RUMi*NDRUH<T6H 

00 1259 J-K 75 »K 76 • 2* 

L»J*23 

REAfi<<0R0H7l C vstr t 1 1 * l*J*L 1 

1759 CONTINUE 

CALL RANGETf E0R|j*2t ISAVEl 

CtaB 

1 

’iso- IEIK2 B.lt, R72> *0 TO 1*0? 

CALL RANSET<K0RuM 2, ISAVEl 
REAOIKORUH?) CVsTRU J »t»K27,K2S) 

•tab 

1 **- ir U 39 .LT.IC 7 *) go To 1500 

CALL RANSET ( K0RMM2 »MORUH < 77 » I 
REA0IK0RUH7) ( VrTK l I » , I »K2* * K30) 

* 

landmark table 

’.Aft* IE* K32.LT.K3 1 1 GO To 1550 . 

’ CALL RANSET I K0rUM2 # KORUN < 75 1 J 
REA0IKDRUM21 ( V 5 TR « I ) » |>K3t »K32) 


-»4*6 3*«# C TEA** LINK INlXfAt. IfATION 

f i J ' u ' . f 


Table B-8 - Partial U stl ng Subroutl ne TRJSUP 


7 

CO 


091*0 

147 

09 1*0 

1*9 

9912* 

I7| 

0*177 

177 

0?|73 

1*3 

Col 73 

177 

09123 

I7S 

0? t 23 

17* 

5n| ?S 

17* 

5*I2A 

)*P 

??!!i 

1 79 
1 An 

05131 

i«1 

“ ; 1 1 7 

l«7 

C5t 3R 

t**% 

37 

IP7 

9*1 3S 

175 

051 3A 
t- JRO 

IAA 

187 

C?t72 

IRQ 

r*»iR2 

17? 

*•> » R3 

|70 

5*1*17 

I7t 

C.-5 77 

17? 

*9!“S 

17?, 

;?i<.7 

IV.i 

P?l 77 

|7r, 

05150 

178, 

9*151 

177, 

9*:«3 

17*, 

?: iss 

!»?, 

C ?»53 

2*", 

**"153 

201, 

C-IS7 

25?, 

99157 

20?, 

CCI57 

29", 


205' 

95 1 57 

29ft, 

95157 

797, 

CC 1*5 

7«r, 

** 157 

75V, 

55 1 A 1 

7 1 n, 

05 1 A* 

2I1< 

9*1*3 

217, 

?5l4** 

2 1 3' 

5* 1 72 

21 Hi 

05 1 7J 

215, 

9*177 

21 A, 

C*I75 

217, 

00174 

21 A, 


t KONFfX ( 313) 
SUBROUTINE TRJ*OP 


. .... It . <RTVPE 

c««*cno standard fOnmon for 
INTEGER V5fRM| 

0 1 Mf NS I ON [r.M3,?l 

DATA ftj /AHVEW I »AHTRAJCC*6HT0RY 

• «hVEH 2 »AHTr/UEC.4HT0*V 


10/0*1/7* 


c 

c 


CALL RT|Hri-*»I T JII # KVeHH 
ghLunt - CGMROt) 

J • A •> *KVrH • I) ♦ 1 

5sZ*- 3 NSZ» fi ° PR * ' * «0NSOL«I.J> ♦ K0NS0L*7tJ) » 

IF C KWCM ,fo» ft MSI • «S77 

CALL DAUXaCVSTR(HT I » , VSTRImYPI > , VsTR < MSUNP > ,NS 2 » VSTR < HypP I) , 
“TRIMBTA8I »VSTR<MKTAe» *VSTR|MWTAB» I 


IP 

21 

, VSTRIMCWFJIOI 

Y5TR<McSEI) » 


VSTR(MVTAn) f VSTRI 
CALL TRIGA( VSTR > 

IF < IC« *En« II ) GO TO 
fF I ICO «HF • 3 ) GO TO 
CALL GPOTAI VsTp(Mcv»E J l I 

I VSTrJNJC,) 

GO TO IS 

|A CALL GPOTa( VSTR<MCphJ|> , VSTR I Me WMJ I *3 I 

| VSTRI MJM| ) , VSTPfMcSMJJ I 

|A I F I M*COM»CQ,0«0p»KVFm»f0*2* GO TO 20 

CALL MCtN|T(V5TRl,MMC0N» ) * VSTR I MMCC« I I » VSTR IhMCaGA ) , VSTRI MMCASA I . 

| VSfRCMVTAnll 

?»* CONTINUE 

IF I NbRNIkvCH) .EG. 9 ) CO TO *CC 


* VST"' HCHCC I > t 
i V5TRIMCRMCI) t 


• Nr I r 

• •-I 

• NE A' 

• NE* 

• WE* 

• •-3 


S 

c 


Call phasea * vsTrcmttpgii , 

| VsTH ( HC N | c I > * 

7 VSTRINALGMI) . 

. J VsTR(MCMALl> t 

*t»? CALL TRAJ I VSTRCMTTPGl) » 

• VStRIMTPI) . 

7 VSTR(MSEV) * tVSTRlMVTAR) ) 

«IR;TE special EVENTS oh Drum 


VSTMIMaNXI | 
V*l TR ( M J Grt I ) 
VSTRlMlNTI ) 
VSTRI MY I I 
VST* I MVTRG) 
VSTRlMYPPI) 


VSTrIMCMLPII t 
VSTRIMLPSII t 
VSTRIMLNYII I 
VSTRlMYPlI I 
VSTRl MY I I * ) 

vstrimoifii , 


• NEW 

• •-I 


3*9 


IF I 
IF I 
N I - 
N2 » 
CALL 
RRITf 


nvemseikveh) »EQ. 0 » GO TO 390 

iscJi 1 * #LT * ? ,AN °' 

M5EVI * ( nvehsf(kveh) • M I - l 
RaNLAS I KoRUMJ? , NDRUM I < VFH-»56 > 
< X0RUM2 J I VSTRI I I i I i HI l 


RTVPE *E0 • HRTYPeO) I GO TO 300 


CONTINUE 
CGPR I I I I ■ rMLUnT 
CALL RTI«E(2»ITjl|,lfVEHI » 
RETURN 
END 


I 

N2 


rno or rONPii ati on * 

T#^$UP SYMBOLIC 

TRJsUP CftOC RELOCATABLE 


NO DIAGNOSTICS. 


07 HAR 7B |**| 151*13 
07 MAR 7B | ** I | 5 I 3 


0 

I 

0 


02** 1 0227 
023Is2)0 
023 1 575*1 


I* 

34 

17 


*1B JCEl.r C 

coelete 




y 


OL MAP.* HOPf, M OPE 


01 Off 


«L *APt* 

I • 

It 

3* 

**• 

5. 

♦ t 

7 . 

fit 

iS: 

Ik 

i**. 

i*. 

I A* 
17* 

I 

!»• 

20. 

21 • 

7? * 
23* . 
?1* 
25. 
74* 
77. 
2R. 

20. 

30. 

3 I * 
32* 

• 33. 

31. 

35* 

36. 

««l SEfi 


Qq 7 SEG 
«Q3 SEjJ 


37. 

31. 

3*. 

10 . 


Table 8-9 


Partial Llstlno Subroutine HOPE 


MO»»f sworf 


o* l>C( f m 


AA 


9 

09 

C 

Cf 

re 

o 

no 

E 

F 

M 

AS 
83 
f 3 
03 
0 


gg 


SCG 

CHN 

SEG 

CHN 

sfg 

sre 


Sfg 

■jrc 

Sr« 

Sffi 

sre, 

SFG 

SCG 

S F« 
SEg 
SFG 
SF6 
SFG 
SFG 
SFG 
SEr, 

SFG 
6LK 
61. K 

CHN 

SEG. 


SfG 

BLK 

CHN 

5EG 


MAIN „ . 

1 

MAIN! -RE TN I T 

2 

INPUT-SrpAR-nAtlti-OOCOl-OECOOE-^lfir-TlHEX^jErCOet 

-OOSL *F-. i QOSPA I r Qf»5nA2.Q0SOA3» iNrCHKtCVpRT , GBPRT t RUNPRT t • 

,!CNP«T ,C0NPPT,A,6.C|0 »E»F tH.BB.CC *00 »SOSP»Tt 

0* Tnp'tS^onST-J « oUMpRC .POTPRC .COVPRC.E0TpKC,TIHPRC.A3»B3if 

| N f Irlf F«-Cnwf /mIs" Jl.PS-SC3D0P 

^pp^m-ac f X YZ 

TP.IP1C-SFTC Oo-SETT«r,-ASS I GN 

SETTAS-. ( ATAMAT »FE ) . . 

GFTnNI>-r,eTSCi.-GETLHl -GET VAU-gCTCAT 

srtnrg-. < oossrt .merge * 

PAPWT-APPKT-SDPRT 

T^SP*?T 

OfiPPT-. csBPRTtlNOPRT.STRpRT.fiODPRTI 
rRDPMC-FN05YM-«< CPOPRl »CROMRG .CRO I HG| 

SFNPRC-r,nnRS-ONOPRC 
STRPPC-siGPRc-OOObS 
« : APRC-I hOPRC 
MCtjpPC — OmPRC 

PATRC-* (BAT I'M .CKALGN.BAFTLL , p A WRT t CKBURN I 

OAT IN - 

00 C01 


maini-super-. ( ; 

FA !R0 ,F* r SUP, POSTROt 

.PSTSUp.rOVRO tC0VSUP.0cR0.6G I . . 

SUOSUP-oCSTOR-OCsur-«<OReOMP.ONBOROI 
DCS TOR 


MAJN3-SUPTRJ-. ITRAJRO.OOI 

OR^^-BrtOT-HATPLV-OUTER-BOOATT-VMASS-SHVENT-MASACC-SHVNTl-l 

•IJACHI A.COESA.USSR) 
solrao-igscon-lopcon 

jgsbRn— lopbrn-s JMULT m i • 

.... - — ■ - * ■ % 

*”oq SEG TRJSUP-TRaJ-TrJOUT.PMaSE-TRIGeR-ENOSTP-OAUX-GPOT-PIMOOI 

-RFADTP-TL^UP; 

-*(Q 01*002.00 3 1 • - ■ • 

QQI SEG ORA6-BODV-mATP L y-OUTER-«OOaTT-VmASS-ShVEMT-MASACC-SHVNTI-» 


•NE» 
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Table C-9 - Partial Listing Subroutine HO^E {Continued) 


•I JACHl A ,t 0 ^SA iUSSM i 

SFG SoL»AO-r«SCON"»i oPCOH 

Sf. 6 Ir, 5 CHN-LOP«ftN-.SJMOLT«t 

scan 

CMM s 

SEG »AT«^- 5 t»P**WO- • CPftOPKO.LI 

sf", ps(r>sur’«iT« jppn-sTTAf'E-'iv*;.* (vrcrnr*act ate ,mm» 

<EG r OV*-**tTr*AJ-*<PPi.PCM‘f»OTAT.MABAT,P*INtT> 
fttK STt/Pf; 

C M N ft 

•vfg m a t N***si*pnn»« t r»Of>p> , a a7» 

v*r, f>r.*-,up-fl(, 7 -r f.R*jp»r .r AtNt r ♦ ,f )PV*Ta,ooy*«tE 7 »CCT» 

s-r, f a;« t ■»*-v'»TM-rrf'Lp ff- : nt««p i -root 
SET, ftUMC AI.-orCAr-*T$r-* < r »ANf,f ♦00rL€*,*AN*tAT»AN6ur»*nA*T»fttrilAC»t 
FPMACC} 

Sr^ 0.nS|t**Tp?»fiCTT»* 1 p»TAI.N»| NPA»T t e*PA0o ,AU.O#*OCCUl,T ,t 

vnisr , v.ir**Mr,,srx mt , pi»nA<f .tM*Af>R » Tapped > 

Srr, pRT J AL-* t C,V:?tT i av«T . COMPT ) 

SfC, nfSTfiH- r nVFC-lMf’t*LT-PW;S 

si.*' OCST^M 


lJ?A 0034 


C‘<* * 

5f -r h a I fiv-SiiPc » r nnwn , t T t 
S!<5 C»DS'JP*C RDtte>j-C«OSAV 
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pROcGSSro MAP 
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o?*oa**7 

otsm%2 

073777** 

9?3/73*| 

07*073*7 

« 7*01*37 
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07*0**37 

97*9*773 
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LjNK 
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